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Zur Berechnung von Gitterströmungen bei einigermaßen 
großem Schaufelabstand 


Herrn R. v. Mises zum 70. Geburtstag 
Von A. Betz in Göttingen 


Zur Berechnung von Gitterströmungen verwendet man vielfach das durch Diagramme gegebene Feld 
einer Wirbelreihe, bei der ein Wirbel fehlt. Dieses Feld wird durch eine Reihenentwicklung dargestellt, 
von der man bei großem Schaufelabstand nur das erste Glied braucht. Weiterhin wird gezeigt, wie man 
aus diesem Feld einer Wirbelreihe, das Feld einer entsprechenden Reihe von Tragflügeln mit gegebener 
Zirkulationsverteilung ermitteln kann und wie es sich auf die Gestalt der Flügel auswirkt. 


For the calculation of cascade flows tha field, given by diagrams, of a series of whirls, whereby a 
whirl is lacking, is often used. This field is represented by a development of series, of which, at a great 
distance of the blades, you need only the first term. It is further shown, how you can ascertain out of this 
field of @ series of whirls the field of a corresponding series of bearing wings with a given distribution of 
circulation, and what influence this has on the form of the wings. 

A la supputation de courants Lreillisses on emploie souvent le champ, donne par des diagrammes, d’une 
serie de tourbillonnements, ou un tourbillonnement manque. Ce champ est reprösente par un developpement 
mangue dont, & une grande distance des aubes, on n’a besoin que du premier terme. Ds plus, il est montre£, 
comment de cechamp d’une serie de tourbillonnements on peut trouver le champ d’une serie correspondante 
d’ailes portantes avec une distribution donnee de circulation, et quelle injluence il a sur la forme des ailes. 


aa BbIyucaeHun pemmeToYHBIX TeyeHHHl HePeAKO VHOTPeÖNAeTCcA 3a 1aHHOE AUATPAMMAMH NONE 
BUXpeBOü HOCHeNOBATENIBHOCTH C ONHHM OTCYTCTBYIINHM BUXpeM. Il 91e 3TO mpencragınerch 
PaA3103KeHHeM B Pal, OT KOTOPOTO Tredyerca, IIPH 6015IHOM PACCTOAHHH JIOHATOK, JIHMB Ile pBbrä 
yıleH. Jlaree HOKA3aH0, KAK U3 3TOTO IOAIA BUXpeBOHÄ MOCHENOBATENIBHOCTH MOIKET ÖBITB HOIYYCHO 
HOAIe COOTBETCTKBYIINETO PANAA KPbLTbeB C AAHHBIM PAcıpenIeHHeM IHPKYJIANUM MU KaK OHO 
HPOABIIAeTCH HA V0OPME KPbIJIbeB. 


Bei weit stehenden Gitterschaufeln kann man für die Berechnung der Strömung von 
einer alleinstehenden Schaufel ausgehen und den Einfluß der Nachbarflügel als Korrektur 
einführen. Dieser Einfluß rührt hauptsächlich von der Zirku- 
lation um diese Flügel her. Als Grundlage für die Berechnung y 
dieser Störung durch die Nachbarflügel dient daher das Strö- 
mungsfeld einer Wirbelreihe, in der ein Wirbel fehlt (Bild 1). 

Dieses Strömungsfeld ist berechnet und in Form von Schau- 
bildern dargestellt). Für manche Zwecke, insbesondere, wenn 


man die Verteilung der Zirkulation längs der Schaufeln be- hr 
rücksichtigen und deshalb die Störungsfunktion über die r 
Schaufelerstreckung integrieren muß, wäre es bequemer, ein- 

fache formelmäßige Darstellungen dieses Störungsieldes ver- Y 

wenden zu können. —>XxX 


So lange sich nun das Gebiet, in dem die Störung dar- 
zustellen ist, nicht erheblich über eine halbe Gitterteilung von 
der Stelle des fehlenden Wirbels entfernt, kann man die 
Störungsfunktionen durch rasch konvergierende Reihen dar- 
stellen, von denen meist das erste Glied allein ausreicht. 

Legen wir in einer komplexen z-Ebene (= x +iy=reiv, 
Bild 1) eine vollständige Wirbelreihe mit der Teilung £so, daß 
die Reihe in die imaginäre Achse und ein Wirbel in den Null- 
punkt fällt, so geht durch die konforme Abbildung & = e’**lt 
ein Streifen der z-Ebene von der Breite der Teilung t in die a 
volle &-Ebene über?). Ist die Zirkulation um die einzelnen 


1) A. Betz, Diagramme zur Berechnung von Flügelreihen, Ing.-Arch. 2 (1931), 5.359. Die gleiche 
Funktion spielt auch in der Theorie des Paramagnetismus eine Rolle, sie wird da als Langevinsche Funktion 
bezeichnet. 

2) A. Betz, Konforme Abbildung, S.166. Springer-Verlag: Berlin/Göttingen/Heidelberg 1948. 
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Wirbel der Reihe 7, so erhalten wir im Punkte ©&=--1 der ö-Ebene einen Wirbel mit der 
Zirkulation Z'und im Nullpunkt, der dem Punkte z= —» entspricht, einen Wirbel mit der 
Zirkulation —//2. Darnach wird das komplexe Potential des Strömungsfeldes 

IE nE — 1 Ib 
2 ni Ye Tor 


8, +1ıY, = In (erle ierat) e 


Das komplexe Potential eines einzelnen Wirbels mit der Zirkulation /’im Nullpunkte der 
z-Ebene ist 


Re 


Der an sich willkürliche Faktor 2r/t von z ist dabei so gewählt, daß für z—0 die Differenz 
der beiden komplexen Potentiale ebenfalls —0 geht. Ziehen wir (2) von (1) ab, so erhalten 
wir als Feld der Wirbelreihe, bei der der Wirbel im Nullpunkt fehlt, 

n2lt __ „rzft . ı 
en - S ar Br RE | 
2700 2nz/t 27 nz/t 


+ = 


Diesen Ausdruck können wir in eine Potenzreihe nach z entwickeln. Es ist nämlich, so lange 
u >27 


Ein u/2 u? ut u® 

In Me Saar HAT nm trPfeat nA ES EA EN 

wobei B,„ die sogenannten Bernoullischen Zahlen bedeuten: 

1 1 1 1 
B= 5: B=—nBe= gg Bm age 
Damit wird 
PER "ik h Kar 1 er rt 2 | 

en 1), RD) 


Für die Geschwindigkeitskomponenten «’ parallel der Richtung die Wirbelreihe und v’ senk- 
recht dazu ergibt sich 


he = er] u 1 ei | 
ur ee 5 + 995 au 2 EEE ERST A): 
Die Trennung in Real- und Imaginärteil macht keine 
y Schwierigkeiten, wenn wir die Polarkoordinaten r und & 
benützen, da ja 
N n 

N =) — ns (cosnp-+isinng) . . . (8) 

0) t i 

* b dl’ ist. 


Ist die Anströmung zum Gitter in großer Ent- 
fernung durch die Geschwindigkeitskomponenten «, und 
IF v parallel und senkrecht zur Gitterrichtung gegeben, so 
hat die Abströmung in großer Entfernung hinter dem 
= X Gitter die Komponenten 4,—=u, + I/t und v. Die Anströ- 
mung auf einen einzelnen Flügel hat nach Wegfall des 

Einflusses der Nachbarflügel die Komponenten 


= (u + W)/2. - a 
und v. 


Geht man von der Strömung um einen Einzelflügel 
von der gleichen Zirkulation ]J' bei dieser Anströmung 
aus, die man mit bekannten Verfahren berechnen kann, 

Aue so ist hieraus in erster Näherung sowohl die Form der 

Flügel yr<&r> als auch die Verteilung der Zirkulation y<@r> 

gegeben®). Aus der Form können wir zu jedem Aufpunkt mit den Koordinaten z=2+1y 
Entfernung r und Richtungswinkel 9 von dem laufenden Punkt z#=2r-+ tyr der Schaufel 


3) <z> soll Funktion von & bedeuten. 
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also z—zr=re'v angeben (Bild 2). Die Elemente dxr der einzelnen Schaufeln haben die 
Zirkulation y<zp> dzp. Demnach ergibt sich der Einfluß dieser Elementarwirbel der Nachbar- 
flügel auf das Potential 
Yw’—-ı® | 
dP — ıdd = = KE— REI NS AUF 0. 3...:.&00 ee (10). 
Dabei ist ® und W’ nach (6) zunächst als Funktion von r und @ gegeben, aber da r und o 
selbst für einen gegebenen Aufpunkt z sich durch x7 ausdrücken lassen, so ist auch ®’ —ı WW’ 
als Funktion von x7 darzustellen. Die Wirkung der ganzen Nachbarflügel wird darnach 


a a er Ba ee (11) 
und entsprechend 
Wr Yy ar PRONS EIG VER ILTI 5. 


Die Integrale erstrecken sich dabei über das ganze Flügelprofil von der Eintrittskante bis zur 
Austrittskante. An sich ist die Darstellung von 7’ usw. als Funktion von xr verhältnismäßig 
umständlich, aber praktisch kann man sich dabei meist mit ziemlich grober Annäherung 
dieser Darstellung begnügen, da ja durch die Integration Abweichungen von der genauen 
Verteilung des Integranden weitgehend ausgeglichen werden. Eine weitere Vereinfachung 
tritt ein, wenn man sich in den Reihendarstellungen (6) und (7) jeweils auf das erste Glied 
beschränken kann. 

Dann kann man nämlich die Integrale auf statische und Trägheitsmomente der Zirku- 
lationsverteilung y zurückführen. Sind &,, %, die Koordinaten des Schwerpunktes der Zirku- 
lationsverteilung, ferner J,—= [y(Yr— Y) dar und J,— [y(yr— 4)’ der die Trägheits- 
momente um die durch den Schwerpunkt gelegten Achsen in &- und %-Richtung und 
PN [rar — %,) (yr — %Ys) dp das entsprechende Zentrifugalmoment, so wird 

[y @- zr) dar = I(@ — z,) Be all. . eo. 7. (13). 
fr (2 — 2? da —J, — St 2 WW yt+ II®— 5) — (y— y)?+ 21:8 8) (y— Y)l - (14). 

An sich kann man auch die entsprechenden Integrale der höheren Glieder von (6) und (7) 
durch höhere Trägheitsmomente ausdrücken. Die entstehenden Ausdrücke verursachen aber 
verhältnismäßig viel Rechenarbeit und die Reihen konvergieren auch nur in einem beschränkten 
Bereich, so daß der Gültigkeitsbereich nur etwa doppelt so groß ist wie bei Beschränkung auf 
das erste Glied. Deshalb dürfte sich eine solche Entwicklung im allgemeinen nicht lohnen. 
Wenn die Beschränkung auf das erste Glied nicht ausreicht, wird meist die Auswertung der 
Integrale unter Verwendung der graphischen Darstellung der Strömungsfunktion einfacher 
sein. Immerhin kann auch diese Entwicklung manchmal nützlich sein. Der Vollständigkeit 
halber sei sie deshalb kurz angedeutet, wobei wir uns auf die komplexe Darstellung beschränken, 
ohne die verhältnismäßig umständliche Trennung von Real- und Imaginärteil ausführlich 
durchzuführen ®). 

Wir führen die dimensionslosen allgemeineren Momente der Wirbelbelegung 


1/a\R 
M, =7 (2) jer- ZN 3 3 el. bainı ed, 02: (15) 

ein, wobei M,=0 ist. Wenn wir zur Abkürzung 
N N... (16) 


setzen, erhalten wir damit 
he 1 BEA «di | R 10: 1 | 
Ba A eat ar ey £ ae Deabtgf,/. BEHUEEEE VORE SITE WER, | PER 
ie 52 o7tM tt 
oe zur 
+ (06 4524 3002 


1 Y 1 Y, .e.* 
u + | 
v—i0- [w+ ind, [@+in)az Se (18). 


4) Für die Berechnung danke ich den Herren Dipl.-Phys. E. Becker und cand rer, nat. K. Kirde, 
9%* 
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Die Parallelströmung, in der sich der Einzelflügel befindet, hat die Geschwindigkeits- 
komponenten % (9) und v und daher in einem Punkte x,, Y, Potential und Stromfunktion 

DV FuYy, WON UI 

Durch Überlagerung der Störung durch die Nachbarflügel gehen die Werte in diesem Punkte x,, 


Y, über in +9, Y4WV, u+W, v-+%. Punkte gleichen Potentials und gleicher Strom- 
funktion verschieben sich hierbei um die Strecken 


Ben aa a ee jr, 
u? + vu? 

Diese Änderung der einzelnen Punkte stellt eine konforme Abbildung der Parallelströmung 
in die von den Nachbarflügeln gestörte Strömung dar. Wenden wir diese konforme Abbildung 
auf die Strömung um den Einzelflügel an, indem wir auch hier alle Punkte um die kleinen 
Strecken nach (20) verschieben, so erhalten wir die gesuchte Strömung durch das Gitter. 
Waren in einem Punkte x&,%, der Umgebung des Einzelflügels Potential und Stromfunktion ® 
und y, so ergeben sich unter dem Einfluß der Nachbarflügel die gleichen Werte in dem ent- 
sprechenden Punkte 2=2,+%, y=Yy,+y’. Voraussetzung ist hierbei, daß die Verschie- 
bungen x’ und y’ nicht allzu groß sind, was meistens der Fall ist. Durch die Verschiebungen # 
und y’ nimmt auch der Flügel eine etwas andere Form an. Deshalb müßte man genau genommen, 
die Rechnung mit der geänderten Form Yr<&r> und Zirkulationsverteilung y<zr> wieder- 
holen. Im allgemeinen sind aber die Formänderungen nur klein und die Wiederholung bringt 
nur so geringe Unterschiede, daß sich diese Iteration erübrigt. 


Fehlerabschätzungen zum Iterationsverfahren bei linearen 
und nichtlinearen Randwertaufgaben 
Herrn R. von Mises zum 70. Geburtstag gewidmet 


Von L. Collatz in Hamburg 


Funktionalanalytische Methoden beginnen sich in steigendem Maße in der praktischen 
Analysis zu bewähren. Freilich bedürfen sie oft einer Umgestaltung, und es bleibt häufig so 
viel Freiheit in der Art der Anwendung. daß die für numerische Zwecke geeignete Auswahl 
noch einiger Forschung und Erfahrung bedarf. Hier soll an dem Beispiel der Iterations- 
verfahren bei Randwertaufgaben gezeigt werden, wie bei Wahl verschiedener Normbegriffe 
die funktionalanalytischen Methoden zu unterschiedlichen Abschätzungen führen, so daß 
die Frage nach zweckmäßiger Auswahl unter den verschiedenen Möglichkeiten auch hier 
sinnvoll erscheint. 


1. Die allgemeine Abschätzung 
Wir legen einen metrischen Raum ®R zugrunde; die Punkte f, fi, fa, ... von N dürfen 
irgendwelche Elemente, Vektoren, Funktionen usw. sein. Für je zwei Elemente f, /, von R 
soll ein Abstand || —fz|| definiert sein mit den Definitheits-Eigenschaften 


IA-FlI20, IA—Al =0 genau dann, wenn i—=ß, 
und der Dreiecksungleichung für beliebige fi, fa» fs: 


IA-RlSIA-Fl + BR - 


In den folgenden Anwendungen handelt es sich stets um Räume, in denen eine Addition mit 
den üblichen Gruppeneigenschaften definiert ist und bezüglich dieser Addition ein Null- 
element © existiert; dann bezeichnen wir || f—©|| = || f|| als Norm von f und definieren an 
Stelle des Abstandes einfacher die Norm. 

Manchmal ist es zweckmäßig, den Raum N noch als vollständig und damit als Banach- 
schen Raum vorauszusetzen; das soll jedoch nicht allgemein verlangt werden. Die Voll- 
ständigkeit besagt: Gibt es in N eine Punktfolge f, mit lim || m —fn|| =0, so existiert zu 
ihr ein Grenzpunkt fin X mit lim || f— || = 0. 

Nn> © 


’ 


Ferner betrachten wir Operatoren 7’ mit Lipschitzkonstanten; der Operator 7 sei 
für die Elemente f eines Bereiches BEN definiert und es gebe eine Konstante K (Lip- 
schitzkonstante) mit 


ITA-TLISKIA-FI für DieEB. 2.0000. 


De 
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Die Operatorgleichung f—= Tf besitze nun eine Lösung « in ®; dann hat für das Iterations- 
verfahren 
Mr Lu, (Dem een. 2), 


im Falle X< 1 Weissinger!) eine Fehlerabschätzung angegeben 
K 
le. — «|| Sg |W%-ı en: 46 9 2 (9) 


ist der Raum sogar ein Banachscher, so läßt sich die Existenz genau einer Lösung im Falle 
K<1 beweisen ?). 

Diese Formel (3) umfaßt sehr viele in der praktischen Analysis bereits bekannte und 
auch neue Fälle; durch Wahl verschiedener Normbegriffe ordnet sich eine Anzahl der bisher 
bekannten Fehlerabschätzungen bei linearen und nichtlinearen Gleichungssystemen hier- 
unter?); ebenso hat man bei Randwertaufgaben verschiedene Möglichkeiten für die Wahl 
einer Norm, wie im folgenden beschrieben wird. 

Folgende Bemerkung über den Bereich ist manchmal ganz nützlich (vgl. das Beispiel 
in Nr. 6). 

Aus 

I, — wl|=|| Tu — Twl| s Klo — W|| 
und allgemein 
| nr — wm Ss EP. — ml; N. 1, 2 VER) 


folgt im Falle KÄ<1 für O<n<m 
| 


r 


Er 2. K 
= 3 Re) l Un — Un—ı | = 1 5 K 


Um — nl). 


m — m || = z (U+ı — %) 


Wählt man sich ein «,E® mit uw,€® und kann man feststellen, daß alle Elemente @ mit 


r 


zu ® gehören, so ist man sicher, daß sämtliche Iterierte «, in ® bleiben. 


2. Der Operator T bei linearen Randwertaufgaben 
Es sei eine lineare Randwertaufgabe mit der Differentialgleichung 
Lupe role) in Bi 2.2.0... 0.0.) 
und den Randbedingungen 
rs le Asse ie RATE 1 Fr 
vorgelegt. L[w] und U,[w] seien lineare homogene Differential-Ausdrücke für eine Funktion 


u, ..., Zm) Ps 9, Yu gegebene Ortsfunktionen, 7’, gegebene Randflächen des Bereiches B. 
Das Iterationsverfahren lautet dann 


L[w. +1] =pWm+t q in B \ 7 
el), Ile 2, aA (79? 
U [un +1] = auf IE J ( ) ( 
Wir setzen voraus, daß die Aufgabe 
Diele re ineB: Dasein. . (8) 


bei beliebigen stetigen r, y„ stets eine eindeutige Lösung 9 besitzt; dann ist also das Iterations- 
verfahren unbeschränkt durchführbar. Öfter werden wir die Existenz einer Greenschen 
Funktion @(z,, &,) fordern, welche die Aufgabe (8) für y„,—=0 in der Form 


en) = | 0a, Er lan Se Be u €.) 


löst. 
Mit Hilfe einer die innomogenen Randbedingungen erfüllenden Funktion % (die auch 


so oft stetig differentiierbar ist, wie es zur Bildung von L[%] gebraucht wird) kann man zu 
einer Aufgabe übergehen, bei der die Randbedingungen homogen sind (dieser Übergang dient 
nur zur Aufstellung der Fehlerabschätzung; er braucht nicht für die numerische Rechnung 
vollzogen zu werden); es sei v=%-+v, „=4-+,; dann ist eine Transformation 7 durch 


Ba een D altne ia (10) 
1) J. Weissinger: Z. angew. Matb. Mech. 31 (1951), S. 245. 


2) J. Weissinger: Zur Theorie und Anwendung des Iterationsverfahrens. Math. Nachrichten 1952. — 
D. Morgenstern: Beiträge zur nichtlinearen Funktionalanalysis. Diss. Technische Universität Berlin 1952. 


= n ‚ Z. angew. Math. Mech. 
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definiert, da %, und damit auch v, beliebig sind; es gilt dann 


%+ı = Tr und TE Fra ME 
ferner bleiben die Fehler beim Übergang von u zu v ungeändert: 
werimetn IT RE 
Bei gegebenem f findet man also g= T/f durch Lösung der Randwertaufgabe 
L[g]= pf HR ins Bad, [el 0 aufal n are 
= =g-+purr Da] ale Bas Erle gerad 


3. Verschiedene Normbegriffe 


Um die allgemeine Fehlerabschätzung anwenden zu können, brauchen wir noch die 
Wahl einer geeigneten Norm und die A Die in der Funktionalanalysis oft 
verwendete Norm für eine Funktion f . 


Sy der je ee (15), 


wobei über den Bereich B zu integrieren ist, ist für numerische Zwecke wenig geeignet; denn 
es kann dann bei Vorgabe eines beliebig kleinen e>0 und eines beliebig großen M>0O noch 
eintreten, daß die Norm || fl|< e ist und trotzdem |f| > M an einer Stelle in B ist. 

Eine einfache, wenn auch manchmal etwas grobe Fehlerabschätzung erhält man durch 
folgenden Normbegriff, der sich aufstellen läßt, wenn die Eigenwertaufgabe 


Lie] =4zen Bir u» U;l2]) = 0 auf Dr BI a dar ea 
eine im Bereich B— % I‘, positive Eigenfunktion Z besitzt, etwa mit dem Eigenwert A,, 


u 
und wenn eine in B nichtnegative Greensche Funktion @(%,, &;) existiert; dann kann man 
die Norm verwenden 


20 
Z 


If an von . (17). 


Es gilt dann also 


LEST AN 


dieser Begriff erfüllt die an eine Norm gestellten Bedingungen, wobei wir zulassen, daß die 
Norm auch unendlich sein darf. Jetzt ist nach (9): 


Ti = Th je@. AED -REFIE (N) 
IA TA<Ih- All: Iplmaz JE (a )ZE)AE (20). 


Nach (16) ist das letzte Integral gerade —— As ; mithin 


4, 
JE el maxın 
ET 


Diese Ungleichung gilt auch für die obere Grenze des ee Ausdruckes; damit folgt 


| I ZA- Thal ERlA—Ale Sn ee ee 
mit dem Wert der Lipschitzkonstanten 


\aimasın 5Y; 
m 


Die Bedingung Ä< 11äßt sich in diesem Falle bequem prüfen und damit bei Ä<1 eine sehr 
leicht durchführbare Fehlerabschätzung gewinnen; angenehm ist bei dieser Abschätzung, daß 
der Fehler nicht im ganzen Grundgebiet B durch eine Konstante, sondern durch die Form 
der Eigenfunktion Z abgeschätzt wird. 

Der letztgenannte Vorteil geht bei dem Normbegriff 


WI=MaRlf| one nenne... (3) 
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verloren; bei (23) ergibt sich in B eine konstante Fehlerschranke. Wieder sei in B eine 
Greensche Funktion @(z#;, &;) vorhanden, von der jetzt aber nicht vorausgesetzt werden 
muß, daß sie nichtnegativ ist. Aus (19) folgt hier 


IA TAlSK 


All 
K=Max [| ö)p@ldi 2... 22.0 (24) 


mit 


Man könnte daran denken, für bestimmte Aufgabenklassen, bei welchen @ nichtnegativ ist, 
[@(z; &) d&, auszurechnen und |p|nax aus dem Integral in (24) herauszuziehen; man erhält 
B 


dann eine handliche Formel für K; aber eine Vergröberung bei K kann, wenn K der Zahl 1 
zu nahe kommt oder 1 gar überschreitet, die Abschätzung (3) unbrauchbar machen. 
Es lassen sich noch verschiedene andere Normbegriffe aufstellen, z. B.: 


e a / s n on 
Mar] 


oder 


|I/||= Obere Grenze von 
in B 


wobei W eine fest gewählte in B positive, evtl. auch nichtnegative Funktion bedeutet. Bei 
dem letztgenannten Normbegriff wird nach (19) 


Ih Tel hl e@ PREWE > (26), 
also gilt Abschätzung (1) mit 
ar Je@ DIS LACK I 
Te W(«;) 
In Nr. 4 wird hierzu ein Beispiel gegeben. 


4. Beispiel 
Vorgelegt sei 
} Zus ya elzu$I)>0 an renere en d28); 


wläßt sich deuten als Biegemoment bei einem auf Knickbiegung beanspruchten Stab. Striche 
bedeuten Ableitungen nach &. Einige Iterationsschritte ergeben?): 


2. 30230% 
32u = 30 — 31x + 28 
nel) ya 


Der Normbegriff (17) ist leicht anwendbar und erfordert keine Rechenarbeit, da bei ihm die 
Greensche Funktion nicht explizit benötigt wird, liefert aber dafür etwas grobe Schranken; 
mit 


2 
Z= 055. %, Ay =T; Belt, ep} 


rgibt (22) den nahe an 1 gelegenen Wert ER ; damit ln == Be 4,29; mit 
we gereg ut 1—K m—8 
|| %, — us|| = Obere. Grenze von AT% | 0,00555 
er a 7 
055% 


erhält man nach (3) und (18) die Fehlerabschätzung 
7 


| — u| < 0,0238 cos 5 


me a 1. 


3) Ausführlicher, aber ohne Fehlerabschätzung, in L. Collatz, Numerische Behandlung von Diffe- 
rentialgleichungen. Berlin/Göttinger/Heidelberg 1951, S. 179, 
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Genauere -Abschätzungen erhält man, wenn man die Greensche Funktion explizit heranzieht. 
Hier lautet sie 


Ge: 2 ee 
G(z, &) = D Eee (29). 
nn für Eeza<i 
Für den Normbegriff || f|| = Max |f| ergibt Gleichung (24): 
1 
K— Max 1 G(z,E)(1 + &)dE = Max = (7 — 6 2? — #) a e 
ası 1 jzlsı 12 12 
8 ’g ri 
Dieser Wert von Ä ist viel günstiger als der obige Wert = ; jetzt wird =- ICH 
(man erhält bereits für «, an der Stelle &—= 0 engere Schranken als oben für %,). 
221 
Ti 1 — Ug|| = Max |u, — U,| = —— = 0,002 h.(@ 
Mit ||, — wz|| Nas u, — u,| 30640 0,00274 wird nach (3) 
| en 0,00384 
1% — use” H O4. 


Noch besser wird die Abschätzung mit der Norm (25), wobei man als Win diesem Falle der 
bequemeren Integration wegen ein Polynom wählt, hier zweckmäßig W=1-—-a?. In (27) 
ist das Integral auszuwerten: 


1 


[6«. HABA AE-= 2,(14 ne 


1— 22) (14 — a2? — #), 


J 30 
‚ 3 14——ı 7 4 RR} RK 7 
es wird Ä Fe 30 — j5 und man erhält den günstigen Wat u re 
Mit = ibt sicl it 
it | — %|| — 5559 ergibt sich somi 
27 # 
e u ea 389 (1 
I, — u| < 5o1en x?) = 0,00382 (1 — =?) 


Eine entsprechende Rechnung mit der Funktion were x? (mit K=3)) führt zu der 
Abschätzung: = 


1 
u — u| <0,00325 1 nn »). 


Diese Fehlerschranke übersteigt z.B. bei —=0 den tatsächlichen Fehler 0,00271 um weniger 
als 20 %. 


5. Eine Klasse nichtlinearer Aufgaben 
Nun sei an Stelle der linearen Differentialgleichung (5) die Gleichung 
Liv] FPlzsu)=V0 in Be 0000 20er 


vorgelegt, während die Randbedingungen dieselben wie in (6), also linear sein mögen. L[w] 
habe die gleiche Bedeutung wie in Nr. 2, sei also wieder ein linearer Differentialausdruck; 
F (x;, u) sei eine gegebene Funktion der &,,...., %,, % mit stetiger partieller Ableitung nach u. 
Die Iterationsvorschrift lautet hier 


L[u.+1] =. F(x,, u.) —=() ineB | (n 
U [u41]=y, auf 24 
Wieder sei die Randwertaufgabe (8) bei beliebigen stetigen r, y„ eindeutig lösbar und es 
existiere eine Greensche Funktion, welche die Aufgabe (8) für y„,—=0 in der Form (9) löst. 
Wie in Nr. 2 wird mit Hilfe einer die inhomogenen Randbedingungen erfüllenden Funktion % 
gesetzt uv=U-+v, W=U+4%. Demgemäß ist die Transformation g=7f definiert durch 
L[g] + F(x,u+f) + [Lul=0 in B, T,lg))=0 al, ... (32) 
Es gilt dann %,4ı = 7, v= Tv. 


— 0, 1, er) ea 
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Zur Festlegung der Lipschitzkonstanten X bilden wir mit 9, = Tf, % = Tf: 


Lg — 9] =F(x,% +) —F(o, u+f), U,ln—gl=0 :- - . - (33). 
In einem bezüglich u konvexen Bereich $ des z,- :-:-2,-u-Raumes, der eine Lösung % und 
die w„ enthält, möge F einer Lipschitzbedingung genügen 
Fa, 2’ —F(2,2*)| <N(®)|2e—2* :- - - : 20.0. .(34); 
dann ist also mt A=f —f, H=Tfı — Tf 
BEE eNmenepr UTdBy 0 aut l!, .. .. ‚Rlnar(35). 


Nun können zwei Wege begangen werden. 


1. Weg. Es wird die Greensche Funktion benutzt; es sei W(«;) wie in Nr. 3 eine in B 
positive (eventuell auch nichtnegative) Funktion, und als Norm wird wie in (25): 


Pe 2: Ce. 00) 


| /|| = Obere Grenze von 
in B 


benutzt. Wir schreiben (35) in der Gestalt 

L[Hj=#-N.|r| mit Jel<ı, 
und nach (9) 

H = /#(6&) G(2;, &) N(&;)| k(&;)| dE;; 
daraus: 


IH < /I@(o;. &) N (&) h(E;)| dE, < all fe; E)NE) WE) - - - - BMI; 


als Lipschit zkonstante ist daher 
j [ea &) N(&)| WE) dE; 


K — Obere Grenze von ? ee. 7:3 (06 
in B W(«;) ( 


verwendbar. Diese Formel enthält (27) als Spezialfall. 


2. Weg. Wir benutzen den Begriff der Randwertaufgabe monotoner Art in einer speziellen 
Fassung®), wobei wir uns weiterhin auf reelle Größen beschränken; wir sagen, bei der Aufgabe 
(30), (6) liege monotone Art vor, wenn aus 

LelzL[ul@inB Je =dU,lu) = sau, rar N 
folgt v<w; dabei sollen v, w zwei Funktionen mit stetigen partiellen Ableitungen so hoher 
Ordnung sein, wie sie bei Z und U, benötigt werden. 


Nun sei 2(%;) eine Funktion mit 


LEI = A207 an BIT, RI 0 ana, a 0 3941 (40). 
Dabei sei & eine Konstante. Unter unseren Voraussetzungen folgt dann aus 
|Z[y]|<D=const in B, eelauT, =... . . (4) 
die Abschätzung 
D [2 
WIS— 22.2 ne nenne: 


denn aus 


re <s—D<L[—y] und LysD<L u | 


. (43) 
zD ( 
0,4 ]= lin =0 
erhält man wegen der monotonen Art 
D 
vlsl2lZ (44) 


hieraus ergibt sich y=0 für D=0; aus L[p]—=0, U,[p]—0 folgt also 9—=0; es muß somit 
nach (40)2>0 sein; man kann also in (44) bei z die Absolutstriche fortlassen und hat damit (42). 
Nun wenden wir (41), (42) auf (35) an: 


IHis 


z|N Al max AN ö (45). 
& 


4) Etwas allgemeiner in L. Collatz: Aufgaben monotoner Art. Arch. Mathematik 3 (1952). 
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Ist W(«;) eine Funktion wie beim 1. Weg, so erhält man mit der Norm (36): 


Naar 
| || hr gi EEE 
also ist als Lipschitz konstante 
N max 
a 2 . (47) 


verwendbar. 


6. Beispiel einer nichtlinearen Aufgabe 
Die Randwertaufgabe laute: 


1 
ee ee ee tk: 


u(x) läßt sich deuten als stationäre Temperaturverteilung in einem Stab(Enden = +1 auf 
konstanter Temperatur v=]1, an die Umgebung der Temperatur u—=0 findet Wärmeabgabe 


1 2 i } 8 : 
proportional + — u? statt). Die Existenz einer Lösung u(z) im Bereich 0 <u<[1 ist physi- 
prop 1 5 phy 


kalisch evident und ergibt sich auch mathematisch nachher unmittelbar. Beginnt man die 
Iteration mit einer die Randbedingungen erfüllenden Funktion 


%W=1+all— 2) 


(— 0,394 


und bestimmt a aus u, (0)—=u,(0), so erhält man a— 173 


; rundet man ab zua=—0,4, 


wählt also 
1 
Wr +2), 
so wird 
ee 1 £ R m2 S 4 6 
= 3000 (1831 + 10385 2 + 130 0° + #9). 


Zur Abschätzung verwenden wir etwa die Norm 


|| /I| = Obere az von . (49). 


1058 


F 3 3 
Im Gebiet O<u<[1 ist = —=1+ R % höchstens 7: also setzen wir N —= Dr Die Greensche 


Funktion ist dieselbe wie im Beispiel Nr.4, stimmt also mit (29) überein. In (38) benötigen wir 


Er [ea su - ma -2 227, 


12 
Nach (38) ist daher 
— 1223 B) i B) 
= Mas fo oma en 
und 
107 
ei = 2 ale 
|| — %|| = Mez x 5000131 134 02 — 4 | 3000° 
107 
In dem Streifen «, + 1800 (1— 2°), der dem Rechteck |x| <1, 0 <u <1angehört, liegen nach (4) 


sämtliche Iterierte %,. 
Wir verwenden daher als Raum R die Gesamtheit der in <— 1,1) stetigen Funktionen f(x) 
mit O<f(x)<1. Dieser Raum ist bei dem Normbegriff || fj| = Ze |\f| und daher auch bei 
1,1> 
dem Normbegriff (49) vollständig?); jede konvergente Flementenfolge ist zugleich eine gegen 
eine stetige Grenzfunktion gleichmäßig konvergente Funktionenfolge, und die Folge der w, 


5) Vgl. z.B. St. Banach: Theorie des Operations Lindaires 1932. Neudruck New York 1949, S. 11. 
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konvergiert daher ebenfalls gleichmäßig in <—1,1> gegen eine stetige Grenzfunktion «(z), 
welche Lösung der Integralgleichung 


1 
y 1 
u(®) [6 SIKAG) Hi u(E)|dE 
—] 
und damit auch der Randwertaufgabe (48) ist. Im Bereich 0 <u<[l ist die Lösung der Rand- 
wertaufgabe zugleich eindeutig; ob die Randwertaufgabe sonst noch Lösungen besitzt, wird 
hierdurch nicht entschieden. Nach (3) gilt dann 


107 > 0,0595 


Im ul 300°; 


oder 
107 i 
1 7 ma z 
ET 
Engere Schranken für die Lösung kann man durch Fortsetzung der Iteration erhalten, be- 
quemer aber durch Anbringen kleiner Änderungen; so kommt man z.B. zu der Näherung: 
un = 0,608 + 0,352 x +0,04 x! 


|, — ul < 


dann wird 
u, = 0,6086400 + 0,3502080 x? + 0,0382507 x! + 0,0027712 28 
+ 0,0001257 x® + 0,0000044 z1°; 
man erhält || «, — w|| = 0,0035478 
und damit |w—w| <0,005913 (1 — 2°). 
Mit der Norm 


f(&) 

— M - 
ee 2 
2) 


| = 


wie im Beispiel Nr.4 erhält man auch hier in der Mitte des Intervalls <— 1,1, günstigere 
Ergebnisse; es wird 
11 K 11 


— —— — 22; FM =VU, > 
K E 2; || — %|| = 0,000765 


und 


2 
I — u] <0,00109 (1 — ). 


E41 


Der Fehler von u, beträgt also an der Stelle &—=0 höchstens 0,28%. Mit der Norm (23) wird 


3 
K —7und man erhält wieder etwas gröbere Schranken, nämlich 


|, — u| < 0,00192. 


7. Allgemeinere nichtlineare gewöhnliche Differentialgleichung 


Ähnliche Betrachtungen lassen sich auch in verschiedenen Fällen durchführen, bei 
denen die Nichtlinearität der Differentialgleichung sich nicht nur auf «%, sondern auch auf 
Ableitungen von « erstreckt. Es seien die gewöhnliche Differentialgleichung k-ter Ordnung: 


eu) onen.» (50) 
und die Randbedingungen an zwei Stellen a, b: (b>a): 
U,[u=y, Be... (bl) 


gegeben. Hierbei sei k>s, L[u] ein linearer homogener Differentialausdruck in u, die Rand- 

bedingungen linear und F eine gegebene Funktion seiner Argumente. F genüge einer Lipschitz- 

bedingung: 

IE u 2, UM) — Elz, u, ur... u) S I Ale) |u? — u | mit So) a. 
=0 52), 

ferner sei die Randwertaufgabe 


L[ul=r(e), Dei] V; N N ) ee (7) 
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bei beliebigem r(x) eindeutig lösbar mit Hilfe der Greenschen Funktion @(z, £) in der Form 


b 
Schließlich sei diese Formel s-mal nach & differenzierbar 
u (x) = fe (2, E)r(E) dE (7=:.0,,1, 5 Bere . . (55) 
mit 
oa G(z,&) 
/(0) ee SM Na 2r 2 
G9 (x, &) Eu 
Die Transformation g= T/f wird dann definiert durch 
La sk en. D.Wl=% k=li...,%). Se . (56). 


Für Funktionen f, m, 9% mit = Tf; (=1, 2) undA=f\, — fi, H=Tf,—- Tf gilt dann 
L[H)=F(a ff AM) - Fo ef. h®), U[LH]=0 


oder 


HEBT A HE Sr BI | 


b 8 | Fa rc (57). 
<Ier HZ, Al Ar (Sl ds 
Nun verwenden wir als Norm einer s-mal stetig differentiierbaren Funktion 
1 Ss 
— —— (0) NET I. 
I\fil = Obere Grenze von gr; > A, (a) |f® (@) | (58), 


wobei W(zx) eine in a<a<<b positive (evtl. auch nichtnegative) Funktion ist; dann wird 


> Aa) HR (AI <IA 2 A. (8) Siena, SW de 


Als Lipschitzkonstante X der Transformation 7 ist daher 


8 b 
K = Obere Grenze von waz+a| \G(z, &)|.W(E)dE . .». . 60) 


ina<e<b 
verwendbar. 


8. Spezialfall einer gewöhnlichen Differentialgleichung 2. Ordnung 
Als Beispiel greifen wir den Fall einer Differentialgleichung zweiter Ordnung: 
Lv] EEE ee ol 
mit den Randbedingungen u(a)—=u,, u(b)=u, heraus. Wählen wir etwa a=0, so lautet die 
Greensche Funktion: 


56-9 für 0<a<t 


ee) für ) 


G(x, &) = (005, 


Der Einfachheit halber wollen wir A, und A, in (52) durch Konstanten ersetzen. Tritt in der 
Differentialgleichung (61) w nicht explizit auf, ist also A,—= 0, so kann man für W(z) die 


Eigenfunktion sin" wählen und hat wegen 
b 
| ade de 
AR, 
ET: mar Cd > a: (63) 
b 
nach (60) die Lipschitzkonstante 
4,6 
K= 2 2 > (64). 
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Ist jedoch A,+0, so würde für diese Wahl von W(x) keine endliche obere Grenze in (60) 


existieren; wir wählen dann etwa W(xz)=1—.az mit z—= au —_ a) und a< 4; dann ist W>0 
inos<sz<sb. 
Mit 
b 
[6 (@. 9 W(e) de (1 Een) 
k,—= Max ? Max = 2 2 
Tozes W(«) ges 2 1— az (65) 
b 6 h . 579] 
je (&, © W(8) ae 2:45? 
k, = Max ° — Max b—- = 
e 0<a<b W(«) 0<z<1/4 —.&7 


hat man dann die Lipschitzkonstante 
A Ra le 3 ee intae na... (66); 
wir geben einige Zahlen an: Für «=0 wird W=1 und die Rechnung einfacher als bei #0; 
es folgt: 
Ab: 
8 


Ab 


Ri mer Miele danaalnunı + (67) 


für a=0 kann man allerdings mit mehr Rechenaufwand oft schärfere Abschätzungen erzielen; 
für 


& -ı (15 E= \ 33 ) » 2,314 folgt®) mit f = ee = 0,3072 
2 
ee... 
Hat man bei der Differentialgleichung (61) die Randbedingungen (etwa mit a—=0) 
u(0) = Us; w(b)—= u, TE EEE (69), 


so tritt an Stelle von (62) die Greensche Funktion: 


I el erze 
A 


Tritt in (61) w nicht auf, ist also A,—=0, so erhält man mit 


W(a)=sin,, die Lipschitzkonstante 


rs... 270). 


Ist jedoch A,=+0, so muß man wie oben andere Funktionen W(x) verwenden; bei W(x)=1 
erhält man 


2 
a Lan. 0 
und?) bei We) =3 —6r +4 2°: 
; K— Aybe : 0,636 4-5 Ayb a rar. (72). 


6) Vgl. F. Lettenmeyer: Über die von einem Punkt ausgehenden Integralkurven einer Differential- 
gleichung zweiter Ordnung. Dtsch.Math.7 (1944), S.56— 74. Lettenmeyer bewies Konvergenz des Iterations- 
verfahrens bei Ausgehen von einer linearen, die Randbedingungen erfüllenden Funktion u,(x) gegen die ein- 

2 


deutig bestimmte Lösung der Randwertaufgabe im Falle En +4,23 b<1. Der Vergleich zeigt, daß die 


- R 
optimale Konstante “ auch hier bei (64) für A, = 0 auftritt, daß aber im Falle A,-+0 in (68) der Faktor von A, 
ungünstiger und der von A, günstiger ist als bei Lettenmeyer. e 

?) F. Lettenmeyer: a. a.O. hat bei den Randbedingungen (69) die Bedingung Een, age A,b<l. 
Die optimale Konstante 3 tritt auch hier für A,=0 in (70) auf; im Falle A,=-0 ist in(72) wieder (wie in (68)) 
der Faktor von A, ungünstiger und der von A, günstiger als bei Lettenmeyer. 
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9. Beispiel 


Es sei 
L[uJl=—- uw’”=F(e, uvwW)=6s2— uu; u) sul) Er 
gegeben. Aus «= 0 erhält man durch Iteration (nach (2) mit (56)) 
1 
W=C— 0, Un — 10 (202 x — 175 — 420° +15 2°); 


die Unterschiede sind 


ai (2 — 2°) (8 — 27 22 + 15 9) 


U U = 


1 
+ = — (8 — 105 (1 — 2?)?). 
: 510 (8 — 105 ©? ( ar) 
Es tritt wieder die Schwierigkeit der Abgrenzung eines konvexen Bereiches auf, der eine 
Lösung und die Näherungsfunktionen enthält. Wir gehen von dem groben Bereich 


er zz urzslelar 


aus. Der Bereich wird nachher schrittweise mehrmals verkleinert, und damit werden die 

Fehlerabschätzungen stark verbessert; aber zunächst verwenden wir als Raum R die Gesamt- 

heit der in <—1,1> stetig differentiierbaren Funktionen f(x) mit |f|<1, |f|s3. Bei dem 

Normbegriff || {|| = Max (4: |f|+B-|f'|) mit 4, Bals positiven Konstanten ist dieser Raum 
<—1,1 


Uy — Us 


vollständig und es sind für jede konvergente Elementenfolge f„ zugleich f„ (x) und f, (x) gleich- 
mäßig gegen stetige Grenzfunktionen f(x) bzw. f(x) konvergent®). 


IEANS az: 
Mit 


vwd, = lu SA, 1 


ou ou’ 
hat man nach (67) (die Rechnung nach (68) bringt im Endergebnis nicht viel mehr, da die 
Lipschitzkonstante nur wenig kleiner wird) 
3 1 7 Kar 
8er 2 STL—K. 
wir arbeiten also zunächst mit der Norm: 


IS = Max @l/I+ PD: 


| oF oF 


för 


Vorläufig schätzen wir grob ab 

|ü — u,| < 0,0075, Im, — %’| < 0,0382 

| u — || < 3:0 0075 + 0,0382 — 0,0607. 
Zugleich sehen wir, daß der oben zugrunde gelegte Raum AR nicht nur “,, 4, %,, sondern nach 
der Abschätzung (4) auch alle weiteren «, und wegen der Vollständigkeit auch eine Grenz- 


funktion v(z) enthält, derart, daß in <0,1> gleichmäßig u, gegen uund u, gegen w’ konvergieren. 
Diese Grenzfunktion «(x) ist Lösung von 


ua) = [62 9P(&, we), w(8)) de 


und damit Lösung der Randwertaufgabe (73). Wieder ist in RW die Lösung der Randwertauf- 
gabe eindeutig; aber ob sonst noch Lösungen der Randwertaufgabe existieren, wird hierdurch 
noch nicht entschieden. 
Nachdem so die Frage nach der Existenz einer Lösung u(x) beantwortet ist, liefert die 
Fehlerabschätzung (3) unmittelbar 
| u, — u|| < 7 - 0,0607 = 0,4249; 
es ist somit j 
3 u, — ul + | u’ — w| < 0,4249, d.h. 
u, — u| < 0,142,  |ü’ — w| <0,425 und damit 
|% | <|w,| + 0.142 < 0,535 
|w| < |us’| + 0,425 Ss 2,463 . 


a z.B. St. Banach: a.a. O, 


( \ 


ae er Föppl, Ein Beispiel zum Mittelwertsatz der ebenen Elastizitätstheorie 12 


=] 


Diese Werte nehmen wir als Ausgangswerte beim 
2. Ansatz: 


A,—= 2,463, A=0535, K- = Et 


der jetzige Wert —1,36 ist schon wesentlich günstiger als beim 1. Ansatz. Mit der 


Norm 


f) 


führen wir dieselbe Rechnung wie beim 1. Ansatz durch, wobei wir nur jetzt etwas sorgfältiger 
die Funktion 


If = Ma 22031] +.0,535 - 


2,463 - |, — %,| + 0,535] u’ — wy | 


berechnen und ihr Maximum bestimmen. Man kommt so zu engeren Schranken für vu und w 
und damit abermals zu kleineren Werten von A, und 4.. 


Beim 
DE AmSaz 
A, 2,110, 4,=0,4097 also K = 0,469, = 0882, 
und der Norm 
Ifll = Max (2,11 |f| + 0,409 : ||) 
0<a<s1ı 
erhält man 
2,11: |u, — w| + 0,409 - | ug’ — w| < 0,0210 , 


also 
lu, — u| < 0,01; |’ — w| < 0,052. 


Wir brechen hier ab, da eine Fortführung die Schranken nur noch wenig herabdrückt. 


Ein Beispiel zum Mittelwertsatz der ebenen Elastizitätstheorie 
Herrn R. v. Mises zum 70. Geburtstag gewidmet 
Von Ludwig Föppl in München 


Der Verfasser hat in der Sitzung der Bayr. Akademie der Wissenschaften vom 7. Nov. 
1952 über einen Mittelwertsatz der ebenen Elastizitätstheorie vorgetragen, für den der Beweis 
in den Akademieberichten demnächst erscheinen wird. Der Satz lautet folgendermaßen: 

„Schneidet man aus einer Scheibe, die einem ebenen Spannungszustand unterworfen 
ist, wobei Einzelkräfte nicht im Innern der Scheibe, sondern nur am Rande angreifen dürfen, 
ein Loch aus, dessen Abmessungen gegenüber den Abständen des Lochrandes von den Scheiben- 
rändern klein sein sollen, so ist der Mittelwert der Randspannungen längs des Loches, gleich 
dem Mittelwert der Spannungssumme q—=0,-+-0, längs der Kontur des Loches in der unge- 
lochten Scheibe.‘ 

Für eine sehr ausgedehnte Scheibe, die unter gleichmäßigem Zug oder Druck steht, 
läßt sich der Spannungszustand für ein kreisförmiges und für ein elliptisch geformtes Loch 
rechnerisch angeben und daher auch der oben angegebene Mittelwertsatz rechnerisch nach- 
weisen. Hier soll die Gültigkeit des Mittelwertsatzes für eine Lochform nach Bild la, b mit 
Hilfe der Spannungsoptik nachgewiesen werden. Dabei wurde die Scheibe aus Dekorit von 
lem Dicke, 10,4cm Breite und 12cm Höhe über ihre Breite gleichmäßig gedrückt. Dem 
Isochromatenbild la bzw. 3a entspricht der geringere und dem Isochromatenbild 1b 
bzw. 3b der höhere Druck. Um einen möglichst gleichmäßig verteilten Druck in die 
Scheibe einzuleiten, wurde auf den oberen und unteren Rand der Scheibe je eine Schicht 
Löschpapier gelegt und die von oben mit Hilfe einer Schraubenspindel ausgeübte Druck- 
kraft über eine Traverse eingeleitet. Nach unseren Erfahrungen erhält man durch die 
Zwischenlagen aus Löschpapier eine über die Breite der Scheibe gleichmäßig verteilte Druck- 
kraft. Man hätte die Druckkraft durch Zwischenschalten eines Kraftmessers bestimmen 
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können. Mit der durch Einbau eines Kraftmessers erhöhten Vorrichtung wäre aber die Gefahr 
des Ausknickens in Richtung senkrecht zur Scheibe verbunden gewesen. Um dies zu ver- 
meiden, hätte man mit wesentlich geringeren Druckkräften arbeiten müssen, worunter die 
Genauigkeit der Auswertung gelitten hätte. In unserem Fall ist die Kenntnis der Größe der 


1a 


Bild 1 


Drucklast nicht notwendig, da der Mittelwertsatz statt in Spannungen ebensogut in Isochro- 
matenordnungen ausgedrückt werden kann. Damit entfällt zugleich die Notwendigkeit eines 
Eichversuches. In unserem Beispiel ließ sich aus dem Verlauf der Isoklinen feststellen, daß 
die Hauptspannungslinien für einen Querschnitt der Scheibe, der ungefähr durch die Mitten 
der parallelen Seiten des Loches verläuft (s. Bild 2a, b) senkrecht zu diesem Schnitt stehen, 
so daß der Schnitt selbst eine Hauptspannungslinie der zweiten Schar darstellt. Wenn man 


| al 


2a 2b 
Bild 2 


von den sehr geringen Druckspannungen in Richtung diese Schnittes absieht und nur die 
Normalspannungen senkrecht zum Schnitt berücksichtigt, so daß sie allein für die längs des 
Schnittes auftretenden Isochromaten verantwortlich sein sollen, so erhält man aus den Isochro- 
matenbildern la und Ib die Verteilung der Isochromatenordnungen gemäß den Bildern 2a 
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f) 
TER, 
\ 4 a 


Bild 3 


und 2b, aus denen man durch Planimatrieren der Flächen und Dividieren durch den infoige 
des Loches geschwächten Querschnitt die in Ordnungen ausgedrückte mittlere Druckspannung 
zu 

Ia = 3,04 
und 

9% = 4,55 


erhält. Beide Werte sind wegen der oben erwähnten Vernachlässigung der geringen Druck- 
spannungen in Richtung des Querschnittes von Bild 2a und 2b etwas zu klein. 


[ 
© 


TRZyeR: 
N 


+5 


Bild 4 
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Die Isochromatenordnungen längs des Lochrandes wurden sorgfältig bestimmt und in 
die Bilder 3a und 3b eingetragen. In Bild 4 sind über dem abgewickelten Lochrand s als 
Abszisse die Isochromatenordnungen — 6 als Ordinaten aufgetragen, wobei der Nullpunkt 
s—0 dem Scheitelpunkt des Loches entspricht. Die zu Zugspannungen gehörigen Isochro- 
matenordnungen wurden in 3a und 3b mit positiven und die zu Druckspannungen gehörigen 
mit negativen Vorzeichen versehen. Die Kurve ain Bild 4 entspricht Bild 3a und die Kurve b 
Bild 3b. Durch Ausplanimetrieren beider Kurven und Mittelbildung ergaben sich die Mittel- 
werte: 

%,« = 310 
und 
a =dt 


Nach dem Mittelwertsatz sollten diese Mittelwerte mit den oben angegebenen Werten von 9, 
bzw. gq, übereinstimmen. Beide liegen um 2 bzw. 3%, höher als die Werte q, und ,. Man 
könnte die Übereinstimmung noch verbessern, wenn man bei der Berechnung von g, und % 
die oben vernachlässigten geringen Druckspannungen in Richtung des Querschnittes berück- 
sichtigen würde. Dadurch würden sich die Werte von q, und g, etwas vergrößern. Die Über- 
einstimmung ist aber ohnehin schon hinreichend, um die Gültigkeit des Mittelwertsatzes auch 
für ein beliebig aus der Scheibe ausgeschnittenes Loch nachzuweisen, wie es in unserem Bei- 
spiel geschehen ist. 


Einige Probleme Mendelscher Genetik 


R. v. Mises zum %0. Geburtstag gewidmet 
Von Hilda Gevringer in Norton (Mass.) 


Durch Mendels Vererbungslehre wurde in die moderne Genetik die Wahrscheinlichkeits- 
theorie eingeführt. Es ist bemerkenswert, daß auch in diesem Anwendungsgebiet der Wahr- 
scheinlichkeitsrechnung nur eine Häufigkeitstheorie den biologischen richtig Sinn der Aussagen 
liefert. Während für den Mathematiker heute Wahrscheinlichkeitstheorie wohl im weiten 
Maße eine Häufigkeitstheorie ist!), sind wir in den Anwendungsgebieten von einer einheitlichen 
und klaren Auffassung weit entfernt. Gleich wie in Laplaces Essay Philosophique die Wahr- 
scheinlichkeit als der berühmte Quotient definiert wird, während im gleichen Werk in den 
Anwendungen eben dieser Wahrscheinlichkeit von Häufigkeiten und wirklichem Ablauf des 
Geschehens bei Wiederholung die Rede ist, so steht es noch heute in vielen biologischen Dar- 
stellungen: Die Ausgangs-Wahrscheinlichkeiten, die inMendels ursprüngliche Gesetze ein- 
gehen, sind ‚reine‘ oder ‚‚a priori“ Wahrscheinlichkeiten, d.h. Laplacesche Quotienten; 
doch im Fortgang der Untersuchung gewinnt die Häufigkeitsauffassung naturgemäß und 
notwendig die Oberhand: Die Beobachtungsobjekte werden geordnet und gezählt, usw. 

Abgesehen von dieser Inkonsistenz gibt es bereits am Ausgangspunkt, d.h. im Gebiet 
des ‚reinen‘ Wahrscheinlichkeitsbegriffs, Beobachtungen und Vorstellungen, die, gleich dem 
„falschen Würfel“ in ein Laplacesches Schema beim besten Willen nicht passen. Nennen 
wir einige Beispiele: Sowie man über den einfachsten Falleines Mendelschen Merkmals (bei 
dem ‚‚gleichmögliche‘“ Fälle vorhanden sind), hinausgeht, tritt der sog. „crossover-Wert“ 
auf (s. Abn.2), eine Wahrscheinlichkeit (‚„Austausch— Wahrscheinlichkeit“), die nun nicht 
mehr auf Basis gleichmöglicher Fälle verstanden werden kann. Wir werden darüber und über 
die anschließenden Verallgemeinerungen des crossover-Wertes noch ausführlich sprechen. — 
Die Verhältnisse werden noch eindringlicher, wenn man, durch biologische Tatsachen veranlaßt, 
verschiedene Austauschverteilungen für die männlichen und weiblichen Wesen an- 
nimmt. — Schließlich tritt in der Theorie der sog. Polyploide in vielfacher Hinsicht die 
Notwendigkeit auf, über das bloße Abzählen von Kombinationen gleichmöglicher Fälle 
hinauszugehen. 

Einige Arbeiten, die mich durch mehrere Jahre beschäftigten, sind in mehrfacher Weise 
zusammengehörig: Sie gehören dem Gebiet der Populationsgenetik an; sie basieren auf 
Mendels Theorie; die Untersuchung wird durchgehend im Sinne der rationalen Wahrschein- 
lichkeitstheorie geführt, wobei aus einigen Grundverteilungen neue uns interessierende Ver- 
teilungen und Gesetzmäßigkeiten abgeleitet werden; schließlich ist methodologisch jedesmal 
das gleiche Untersuchungsschema befolgt: Nach Definition der Grundverteilungen 
werden Rekursionsbeziehungen für gewisse Verteilungen in aufeinander folgenden Gene- 
rationen gewonnen, diese werden sodannintegriert, und es werden insbesondere die stabilen 


!) Eine Axiomatik wie die von Kolmogoroff ist weder im Gegensatz noch in Übereinstimmung mit 
einer Häufigkeits-Definition und -Theorie. Eigentlich liefert sie die Grundlegung für die allgemeine Mathe- 
matik der Verteilungen. 
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Zustände untersucht, was zur Aufstellung von Grenzwertsätzen führt. Es ist Ziel dieses 
Artikels aus dieser Gruppe von Untersuchungen einiges zu berichten ?). 


1. Die grundlegenden Verteilungen. Wir erklären diese zunächst an dem einfachen 
Fall eines (‚‚diploiden“) Mendelschen Charakters und unter Annahme von Panmixie. (Man 
versteht unter Panmixie [Zufalls-Paarung, random mating] daß alle jeweils möglichen Kreu- 
zungen auch wirklich auftreten und zwar nach Zufallregeln. Wir kommen darauf noch ein- 
gehend zurück.) Mendels erste, sehr kühne Hypothese ist die folgende: Gewisse Charaktere 
— heute Mendelsche genannt — z.B. Blütenfarben (etwa rot und weiß im Falle der Erbsen- 
blüte, rot, weiß und rosa im Falle von mirabilis yalapa) hängen in der einzelnen Pflanze nicht 
nur von einem Faktor sondern von einem Paar von Faktoren ab, den Genen, die wir hier 
mit Rund W bezeichnen. Der Organismus, die Pflanze, insbesonders die befruchtete Keim- 
zelle (Zygote) ist daher vom Typus RR oder RW= WR oder WW. Ein von einer solchen 
Zelle geformter und bei der Fortpflanzung weitergegebener Gamet (Ei oder Sperm) enthält 
nur einen der beiden Faktoren und zwar besteht für jede Wahl die gleiche Wahrscheinlichkeit. 
Ein neuer Zygote der nächsten Generation entsteht dann durch die Zufallspaarung eines 
männlichen und eines weiblichen Gameten, und zwar nach Zufallsregeln, d.h., wie wir sehen 
werden, gemäß Multiplikation unabhängiger Verteilungen?). 

Für jeden Charakter gibt es r Möglichkeiten (Allele), (z.B. die vielen verschiedenen 


' möglichen Augenfarben bei Drosophyla). Wir bezeichnen einen Zygoten-Typus mit (z, %), 


wobei 

ey) ly)z) Bing eyeliähnz,r) el). 
Es gibt daher r(r +1)/2 verschieden Zygoten (verschiedene Genotypen — der Begriff des 
Phänotypus tritt in unserer Untersuchung überhaupt nicht auf). Wir nehmen in dieser ganzen 
Untersuchungan, daß es Sinn hat von diskreten, nicht übergreifenden Generationen 
zu sprechen, und bezeichnen mit w®) (z,y)=w")(y,x) die Verteilung der Genotypenin 
der nte" Generation, dieselbe für beide Geschlechter, wobei natürlich 

UFER. 


5 
Bee Zee nn nee (122), 
= y 
Unter dem ersten Argument in (z, y) oder w(z, y) verstehn wir das von der Mutter, unter dem 
zweiten das vom Vater ererbte Gen. Weiter besteht nach Mendels eben erklärter Annahme 
gleiche Wahrscheinlichkeit für die Segregation des einen oder andern dieser beiden Faktoren 
bei der Gametenbildung. Die Segregations-Verteilung besteht hier also nur aus der 
Alternative mit den Wahrscheinlichkeiten p=q=1/2. Wichtig ist, daß die S. V.nichtvonn 
abhängt und nicht vom speziellen Genotyp. Diese beiden Eigenschaften werden wir 
in viel allgemeinern Situationen stets wiederfinden. Aus den soweit erklärten Verteilungen, 
der Verteilung der Genotypen (V. Z.) und der Segregations-Verteilung (S. V.) leiten wir die 
Verteilungder Gameten (V.G.)ab, p” (x), d.h. die Wahrscheinlichkeit (W.), daß ein Gamete 
der nten Generation den Faktor & besitzt. Diese ist sehr einfach zu finden. Um den Faktor & 
zu übertragen muß ein Zygote a) ein x besitzen, b) es übertragen. Ein Homozygote (2, x) über- 
trägt = mit W. eins, ein Heterozygot (z, y) mit W. einhalb; daraus folgt, daß 


pm (1) = ww(1,1).1+wm (1,2). 43 + wm" (2,1).43 +. + w(r,1)-$ 
also 


1-2 
pm) Zumlsy) wobei SpWd)=1....... (13) 
Y z 


wegen (1.2) gilt. Schließlich folgt gemäß unserer Voraussetzung der Panmixie die Verteilung 
wr+D(z, y) aus den p® (x): 
wr+d(z, y) = pm (z)pm(y), Beueiea. etw (1.4) 


immer für @=1,...,, y=1,...r. Aus (1.3)und (1.4)erhält man nun ineiner Zeile Hardys 
berühmtes Gesetz von der Konstanz der Verteilung der Typen. Denn 


pR+D (a) = Z ward (a, y)= I pwla)py)= Pa), m=0, 1... - (13) 
y Y 


2) Dies mag gerechtfertigt scheinen, da diese Artikel auf mehrere, amerikanische, nicht mathematische 
Zeitschriften verteilt sind; sie umfassen etwa 170 Druckseiten. Diese Arbeiten des Verfassers werden von nun 
an als (G. (a)) usw. zitiert. — Aus Raummangel werden biologische Arbeiten nur sehr unvollständig zitiert, 
d.h. nur, wenn sie in direkter Beziehung zu den referierten Arbeiten stehen. Es sei auf das Literaturverzeichnis 
in (G.(g)) verwiesen. A z 

3) Wir müssen auf biologisch vollständige und befriedigende Formulierungen leider verzichten. Unsere 
Objekte sind mathematisch begriffliche, schematische Modelle des biologischen Vorgangs. „Doch ziehen wir 
es vor, die biologischen Vorstellungen anzudeuten, und nicht einfach mit „Urnenschemen“ zu arbeiten. 

10* 
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und daraus auch 
untl (x, Ye wer) (x, Y), (nl 2, Ä a) 127. KR HER (1:5). 


Es ist hier durchweg Panmixie (random mating) vorausgesetzt worden. Diese ist, negativ, 
durch das Fehlen jeder Selektion definiert, positiv eben durch die Gleichungen (1.3) and (1.4). 
In (1.3) ist ausgedrückt, daß die Proportion der Gameten, p®(x), durch die der Eltern be- 
stimmt ist, eine Proportion, die sich ihrerseits zwischen Geburt und Fortpflanzungsalter nicht 
geändert hat, und daß alle Genotypen ihrem vorhandenen Zahlenverhältnis gemäß (ohne einen 
vom Genotyp abhängigen Unterschied in der Fruchtbarkeit usw.) an der Zeugung der neuen 
Generation mitwirken. In (1.4) ist das Fehlen von Gattenwahl ausgedrückt. Im Gegensatz 
zu all dem würde die Formulierung einer allgemeineren Fortpflanzung etwa von der Form sein 
we, — N nU4 o%,, wobei wir mit wen die W. des Typus ‚u“ (durch beliebig viele 

r,4 
Gene gekennzeichnet) bezeichnen, mit 2 die W., daß zwei Typen x und A zur Zeugung zu 
sammentreffen (was außer vom Zahlenverhältnis der Typen von den verschiedensten Selek- 


tionen bestimmt sein mag) und mit 0%, die W., daß aus dieser Vereinigung ein Typus entsteht, 
eine Verteilung, die nicht das Produkt zweier S. V.en zu sein braucht und auch noch von n 
abhängen könnte. 


2. Vererbung von m gekoppelten Merkmalen. Betrachten wir simultane Cha- 
raktere, wie Farbe der Blüte, Länge des Stengels, usw. Für jeden Charakter sind r Allele 
möglich (es ist unnötig m Zahlen r,, ... r„ einzuführen, da man im Falle verschiedener 
Allele-Anzahlen mit r die größte der r, bezeichnen und entsprechende W.en gleich Null an- 
nehmen kann). Ein Genotyp ist dann mit (x, y) zu bezeichnen, wobei die ‚‚Vektoren“ z, %, 
beziehungsweise die Komponenten 2, ... &, und %, ».. Y"„haben (u=1,...r, y=1,...r, 
@=1,... m)), und x das mütterliche, y das väterliche Erbe bezeichnet. Wieder setzen wir 
wiein Abn.1(2;Y)=(Yy; x), oder ausführlicher (2,2. &n; Ya. Ya) — (Ye Ymsbsauk an) 
voraus. Hier ist es wesentlich durch den Strichpunkt mütterliches und väterliches Erbe zu 
trennen, d.h. den Genotyp durch zwei m-dimensionale Vektoren zu geben und nicht 
durch m Paare von Ziffern. D.h. es ist etwa (1,2; 5,7) von (1,7; 5,2) zu unterscheiden, 
obgleich beide bez. des ersten Merkmals durch das Paar (1,5), bez. des zweiten durch (2,7) 
gekennzeichnet sind. Es gibt dann N,—=4#:.r”(r®-+1) verschiedene Genotypen, z.B. für 
r=m=2ist N,=10, während es N, = [r (r +1)/2]” Zahlenpaare gibt, für r=m=2, N,—9. 
Diese zweite Auffassung reicht aus, wenn man Mendels ‚‚zweites Gesetz‘ annimmt, dem- 
gemäß die m Genpaare voneinander unabhängig sind. -Dies aber ist eine nach dem 
heutigen Stand der Beobachtungen zu enge Auffassung. Die Beschreibung des Erbvorganges 
ist nun die folgende: Im Fortpflanzungsvorgang wird ein Gamet gebildet, der bez. jedes der m 
Faktoren ein Gen enthält, also etwa für m=5: (2%, %9 Ya, %, %). Dabei ist es, wie bemerkt, 
nicht mehr die geltende Theorie, gleich Mendel anzunehmen, daß jeder dieser 2” Kombi- 
nationen die gleiche W. nämlich 27” zukommt, sondern es ist der besonders von Morgan 
beobachteten, sog. Koppelung (linkage) Rechnung zu tragen. Das Grundphänomen der 
Koppelung ist, daß eine Tendenz besteht, die Weitergabe der von den Eltern übernommenen 
Genverbindungen zu bevorzugen, also ein Einfluß der Großeltern auf den Enkel. Z.B. im 
Fall m=2 wird ein Wesen vom Typ (%, %; %ı, %s) die Gameten (2, %) und (%Y,, %) je mit 
einer W. g=(1—.c)/2 abgeben, während für Abgabe der beiden andern möglichen Gameten 
je die W. p=c/2 besteht. Man nennt c, das gewöhnlich <1 ist, den crossover-Wert (auch 
recombination value, oder Austausch W.). Dies ist die erste der in der Einleitung angedeuteten 
Situationen, wo die Laplaceschen 1/4, 1/4, 1/4, 1/4, die den vier Möglichkeiten der Gameten- 
bildung entsprechen, nicht ausreichen. Das gilt erst recht von der allgemeinen Verteilung, der 
wir uns jetzt zuwenden. 

Man sieht leicht, daß man zwar für m=3 durch Angabe der drei crossover-Werte c,,, 
(93; €13 (5 Cj;) die Segregation voll beschreiben kann, daß aber für m >4 eine komplette S. V., 
die wir hier Koppelungs- oder auch Austausch-Verteilung (A. V.) (linkage distribution) nennen 
wollen, eingeführt werden muß (s. G. (a)). Bezeichnet SdieMengedermZahlen1,2,..m, 
Aeine Untermenge von S (die auch gleich 0 oder S sein kann), 4’=8—A, dann sei (A) 
die W.eines Gameten, dessen mütterliche Gene von 4A, dessen väterliche von 4’ 


stammen, und 
KA)= 14), SUA)=1 061 Kine 2er 


Fügt man noch die Beziehung (1.4) hinzu, wobei x und y jetzt natürlich Vektoren sind, so 
ist die Vererbung schematisch folgendermaßen beschrieben: a) Ein Genotyp ist charakterisiert 
durch zwei m-dimensionale Vektoren x und 4; erstes (zweites) Argument bezeichnet mütter- 
liches (väterliches) Erbe. b) Jeder abgegebene Gamet ist von der Form (2, ... 2m), wobei 
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jedes 2; gleich ©; oder y; ist. c) Die jeweilige Auswahl unter den 2” möglichen Gameten erfolgt 
nach der oben erklärten A.V. d) Der neue Genotyp in der folgenden Generation entsteht 
durch Multiplikation entsprechender Gameten-W.en, gemäß (1.4). — Dabei ist angenommen, 
daß für beide Geschlechter die gleichen Regeln und Verteilungen gelten. Es ist klar, daß in 
unserer A.V. alle beobachteten Spezialfälle, wie ‚Unabhängigkeit‘, alle 1, gleich; ‚‚voll- 
ständige Koppelung‘, /3—=1,=1/2, alle andern gleich Null; und weitere Fälle, die wir noch 
besprechen werden, enthalten sind. 

Ehe wir weitergehen, müssen wir an den Begriff der Marginalverteilung er- 
innern: Aus einer m-dimensionalen (arithmetischen) Verteilung P (2, 2... 2m) leiten wir 
eine Marginalverteilung von niedrigerer Ordnung als m ab, etwa die Marginalver- 
teilung 3. Ordnung 34 (23232,) durch (m—3)-fache Summation I... Plz 2 - - . m): 

Zi 2 
Für 9334 (23232,) schreiben wir kürzer p4(24), oder auch 94 (2), wobei A die Teilmenge (2, 3, 4) 
bezeichnet. 

In derselben Weise ist die Marginalverteilung unserer A.V. zu verstehen: Es bedeutet etwa 
l4(A,)=14(4,), wo4,+4,=4A, A< 8 die Wahrscheinlichkeit, daß die mütterliche Gene, die A,, die 
väterlichen, die A, zugehören, abgegeben werden, während für die übrigen Gene jede Möglichkeit besteht. 
Man mag im übrigen die A.V. als eine m-dimensionale Verteilung mit Argumenten g; auffassen, wobei 
&=1 (oder gleich Null) wenn das ite abgegebene Gen ein mütterliches (väterliches) ist. Es ist nun auch 
klar, daß c;£=cgi eine durch (m—2)-fache Summation aus der A.V. entstebende Marginalverteilung ist, 
d.h. es ist die W. dafür, daß das mütterliche Gen i und das väterliche k abgegeben wird, oder umgekehrt. 

Wir wollen nun eine Rekursion für w®) (2; y) oder, was bei Panmixie auf dasselbe heraus- 
kommt für p®(x) angeben und müssen hier aus Platzmangei die Resultate ohne die, im 
übrigen sehr durchsichtigen, Beweise mitteilen. Zunächst gilt als Analogon von (1.3): 


DI HA (2) Be aa nr: (22), 


dabei ist A eine beliebige Untermenge von S, und x und %y so, daß 
ge amenn ae Au = 2, weunned . .... 2... (22). 
Die entscheidende Rekursionsformel ist: 
pn (2) =. LADE DR ZI (2,5) 


wobei die Summe über alle Teilmengen von S zu erstrecken ist. Beachtet man die Symmetrie 
I(A)=1(4’) so gibt (2.3), etwa für m=3: 


PD (21, 22, 2) = 2 [1(0) pP" (21, 22, 23) + (1)Pı (2) p35 (2223) + (2) P2 (%) ps (21 23) 

+ 1(3) P3(23) PVy (21 22)] (2.3°) 
oder für m = 4 in unmittelbar verständlicher Schreibweise: 

pPtdez...2)=2[0)P" (2...) + Para) ++) 

+ (112) pp 9 + 113) pi PP + I(14) pi? p%)] 2.3”). 
Man sieht, daß es die Werte der A.V. sind, die diese Formeln so durchsichtig machen, indem 
sie die Zusammenfassung zu sinnvollen Gruppen nahelegen. Es ist klar, daß wir bei Durch- 
führung der Rekursion (2.3) auch Rekursionsformeln für die Marginalverteilungen 
brauchen. Z.B. treten in (2.3) auch »%) oder p%) auf. Es ist aber leicht zu sehen, daß dies 
kein neues Problem darstellt, sondern daß etwa die Rekursion für p{%, einfach durch (2.3°) 
geliefert wird. 

Die Formeln (2.3) bilden ein System von Quadratischen Differenzengleichungen mit 
konstanten Koeffizienten. Das Ziel der Lösung oder Integration ist, p'*) als Funktion von n 
und der Ausgangsverteilung p‘ =pauszudrücken. Für den Laplace-Mendelschen Spezialfall 
der Unabhängigkeit der Charaktere hat H. Tietze*) bereits vor 30 Jahren unser Problem 
behandelt und sowohl Rekursionsformeln wie auch Limesverteilung für n>oo angegeben. 
Hingegen stammt die elegante explizite Lösung der Rekursionsformel für diesen Fall von 
F. Bernstein°). Die allgemeine Lösung bei beliebiger A.V. ist nicht einfach: Für m—=2 
haben H. S. Jennings®) und auch R. B. Robbins die explizite Lösung gegeben. Für m—3 
ist die von mir (G. (b)) gegebene Lösung noch sehr einfach. Die allgemeine Lösung ist in G. (a) 
enthalten. 


4) ‚Über das Schicksal gemischter Populationen nach den Mendelschen Vererbungsgesetzen.““ Z. angew. 
Math. Mech. 3 (1923), S. 362—393. 
x 5) „Variations- und Erblichkeitsstatistik‘ in „Handbuch der Vererbungswissenschaft“, Bd. 1, S. 1—96. 
6) (a) The numerical results of diverse systems of breedings withrespect of two pairs of characters, etc. 
J. Geneties 1 (1916), p. 97—154. . 
(b) The numerical relations in the erossing over of the genes etc. J. Genetics 8 (1923), p. 393 —457. 
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Was nun die Stabilität betrifft, sieht man leicht, daß jede der in (2.3) auftretenden 
Marginalverteilungen erster Ordnung gleich der entsprechenden in der Ausgangsgeneration 
ist, ebenso wie in (1.5). 


2% (z,) = p (z,) = p(%), Bl... mM) iR — 1,02. Eee an 
Aber schon für m>2 ist ein Gleichgewichtszustand bei endlichem n nur in dem trivialen 
Fall kompletter Koppelung vorhanden (1($)—1(0)=1/2, alle andern l-Werte gleich Null). 
In diesem Fall ist offenbar 


or) (2, zn) == PN 2 (n =0(, il a) re De 6 (2.5). 


D.h. vollständig gekoppelte Gene verhalten sich wie ein einziges. 
Betrachten wir nun den ‚‚Hauptfall“, daß alle crossover Werte von Null ver- 
schieden sind. Dann besteht der Grenzwertsatz: 
lim p® (2)... 24) =P(A)2%): 0: Deu) E re eo: 
Nn>&© 
Im Gleichgewichtszustand, der erst für n— w eintritt, herrscht Unabhängigkeit. (Für Tietzes 
Fall, wo jedes c;;— #, ist unsere Voraussetzung natürlich erfüllt, und es gilt (2.6). Zwischen 
diesen beiden extremen Fällen, — ‚‚alle c;; gleich Null“, ‚Kein c;; gleich Null“, — liegt die 
Möglichkeit von £<m Gruppen &$,, derart, daß innerhalb jedes 8, komplette Koppelung 
herrscht, d.h. alle betreffenden c;; verschwinden, daß aber mindestens ein von Null ver- 
schiedener Crossover-Wert von jeder Gruppe zu jeder andern führt (dabei schließen wir 
auch den Fall der aus einem einzigen Element bestehenden ‚‚Gruppe“ ein). Es ist dann 
S=S, +9, -+:.-8,. Zunächst gilt dann für jede Gruppe 


p8. = Pp8,, E12 2 RE en, 
und für die Verteilung der m Charaktere: 
lim pm, .=Pg2P 28) ER T DES BEFORE EN ZN 


no 


Innerhalb jeder Gruppe kompletter Koppelung bleibt für alle Generationen 
die ursprüngliche Gametenverteilung erhalten. Die m-dimensionale Vertei- 
lung aller Gene konvergiert für n—oo gegen eine Grenzverteilung, in der 
die {-Gruppen unabhängig verteilt sind. 


Eine triviale Verallgemeinerung würde darin bestehen, für männliche und weibliche 
Organismen verschiedene Ausgangsverteilungen w®(z;y), und v®(x;y) anzunehmen; schon 
nach einer Generation von Panmixie gleicht sich dies aus. Hingegen ist es biologisch 
wichtig und nicht trivial (darauf bezieht sich der zweite Hinweis in der Einleitung) für die 
beiden Geschlechter verschiedener Austauschverteilungen zu betrachten (s. G. (d)), 
nämlich 2, und /, und dann ist auch die Annahme verschiedener Ausgangsverteilungen, w), 
und vo für die Zygoten, p® und g(% für die Gameten nicht trivial. Das Resultat ist, daß dann 
dauernd verschiedene Verteilungen 2”), g'” auftreten; die Rekursionsformeln sind ähnlich wie 
bisher,. doch nehmen r® = $ (p® + g®) und Z,=#(l,+1,) den Platz von'p') und ], ein. 
Diese Verhältnisse sind allerdings für m>4 nicht mehr ganz so einfach. Man erhält jedoch 
für jedes m sehr einfache Grenzwertgesetze: Mit n— oo konvergieren p®) und g'* gegen die- 
selbe Grenzverteilung, diedurch r® (z,,....2„)bestimmt ist. Wenn alle ı; = - (a; +c;;)0, 
so konvergieren beide gegen das Produkt der m Marginalverteilungen N (z;) der einzelnen 
Gene (=1, ... m), beziehungsweise im allgemeinen Fall gegen das Produkt der den oben 
definierten Gruppen 8,, ... ‚S, entsprechenden Ausgangs-Marginalverteilungen der r), 

3. Mendelsche Genetik der Autopolyploide. Die Grundverteilungen. Wir 
haben soweit die Verallgemeinerung von Mendels ursprünglichen Ideen auf den Fall von m 
beliebig gekoppelten Charakteren untersucht. Eine zweite wichtige Verallgemeinerung in ganz 
anderer Richtung bezieht sich auf die sog. ,„polysomic‘ (mehrwertige) Vererbung. Hier 
tritt auch, im Gegensatz zu der bisher ausschließlich betrachteten ‚„Chromosom Segre- 
gation‘ die sog. „Chromatid Segregation“ in den Vordergrund. Wie man leicht zeigen 
kann, kommen bei normalen ‚„‚diploiden‘“ Organismen, gleich den bisher betrachteten, beide 
Arten der Segregation mathematisch auf das gleiche Problem hinaus. 

Wir beginnen mit der Chromosomen Segregation für m=1. Das wesentliche Kennzeichen 
eines Autopolyploids ist, daß ein Gamet bezüglich jedes Faktors nicht nur ein, sondern s>1 
Gene besitzt, und der Zygote 2s solche, von denen je s sein mütterliches, bez. väterliches 
Erbe darstellen; jedes Gen kann r Allele annehmen. Bei der Fortpflanzung gibt jedes Elter 
vonseinen 2s Genen sab und die Auswahl dieser erfolgt nach einem — von nund vom Genotyp 
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des Elter unabhängigen — Zufallsgesetz. Ein Gamete wird also durch s Zahlen (z), ein Zygote 
durch zwei Gruppen von je s Zahlen (x; y) gegeben. Es entspricht den Beobachtungen anzu- 
nehmen, daß beim Gameten wie beim Zygoten den verschiedenen Permutationen der gleichen 
Allele derselbe Typus entspricht. Dementsprechend wird man z.B.für s—6 statt (a, a,a, 0, a, @,) 
r+s—1 
8 
wenn man den Zygoten durch sein mütterliches und väterliches Erbe gibt, wobei (2;y)—=(y;&), 
so $&- R(R+1) Zygoten. Bezeichnen wir nun, in etwas ausführlicherer Schreibweise mit 
uw" (atı...atr; ahı...a’r) wobei u + +%,=Ss, Y-+'--Y,=s die W. eines individuellen 
! s! 
le leo yl 


nun (afa,az) schreiben. Man sieht, daß R = ) Gameten in Betracht kommen und 


Typus, so erscheint gemäß den möglichen Permutationen so ein Typus 


-mal und es ist, mit einleuchtender Abkürzung 
s $ Bu 
ID) een. ....-00 


Analoges gilt für die Gametenverteilung p' (at... ar) und es ist iM Ä A) p”(afı ...a?r) die 
W. daß ein Gamete den Typus af: ...a’r hat und 


Do Wake) — 1 ern: (32). 


9) 
Nun die S.V. Es können s Zahlen aus 2s solchen auf s-(,) Weisen ausgewählt 


werden und es wurden, als einfachste Annahme diese 8 Fälle stets als gleich möglich angesehen, 
insbesondere von Haldane”), an dessen wichtige Arbeit unsere Ausführungen (G. (c)) an- 
knüpfen; dies wird als Zufalls-Chromosom-Segregation bezeichnet. Doch haben Beob- 
achtungen Unstimmigkeiten mit der Theorie gezeigt und aus diesem Grund, sowie um größerer 
Allgemeinheit und — in letzter Linie — Einfachheit willen setzen wir unsere S.V. nicht von 
Anfang an als Gleichverteilung voraus. Studium von Beobachtungsergebnissen führt zu der 
Vermutung, daß die Proportion mütterlicher und väterlicher Gene im Gameten eine Rolle 
spielen mag. Wir führen also „(a1 ein, als die W. dafür, daß eine 
bestimmte Gruppe von & mütterlichen Genen abgegeben wird. Da eine solche 
mit irgend (s—o) väterlichen Genen kombiniert sein mag, so ist 


I>()n=- DI) r=1 ee 05); 


‘=0 a=0 

Für A,=1/S ergibt sich Müllers und Haldanes gleichförmige S.V. Auf diese Weise sind, 
wenn s=2u oder s=2u-+1, u verfügbare Parameter eingeführt. Es sei noch folgendes 
bemerkt: Obgleich wir mit Absicht alle Definitionen und Formeln in möglichster Ähnlichkeit 
zu denen des Abn. 2 gegeben haben, ist doch die Situation hier eine ganz andere als dort: 
Dort gab es r” verschiedene Gameten, hier R=(" 1 r 
polyploide, s>1, m=1 kann jede Gruppe von s Werten abgegeben werden, während beim 
früheren Problem, s=1, m>1, ein ‚‚Vektor‘ abgegeben werden muß, d.h. jede der m Kompo- 
nenten muß einmal vertreten sein. 

4. Fortsetzung: Rekursionsformeln. Grenzwersätze 

Wir suchen wieder eine Rekursionsformel für p®. Es sind die Marginalverteilungen zu 
definieren. Es bezeichne (in irgend einer Generation) p(aa®...a”r) die Marginal-W. der 
Ordnung a, für die Kombination (afı... ar), (4 +: %,—=a, aSs), während die übrigen 
(s—a) Gene beliebig sind. Dies ist z.B. die W. dafür, daß die ersten z,-Plätze mit 
dem Allel a,, die nächsten z, mit a,,... die letzten der a Plätze mit a, besetzt sind. Die 
gleiche W. besteht natürlich für andere a Plätze, wenn die a, und #, dieselben sind; aus 
diesem Grund können Indizes in den Marginalverteilungen weggelassen werden. So ist etwa 
fürs, r=2: 


). In der Segregation der Auto- 


?(4%) = P(4%) +3 Plaia) +3 pl) + p(aı@) 
die W. dafür, daß das erste Gen a,, das zweite a, ist, während die andern drei beliebig sind 
und es gilt (wenn wir dies eine Mal Indizes schreiben) 


Pı2(4ı%) = Pıssald%) + 2 Pızss( ia) + Pızsı(Cı%) = P123(1%) + Pıas (a @2), 
?) „Theoretical Genetics of Autopolyploids“, J. Genetics 22 (1930), p- 359— 372, 
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Unsere allgemeine Rekursionsformel ergibt sich in der Form: 


$ 
Z\[2 z 
Pen lapat. II Zla)la)- (armer. mama... 
a=0 %ı %, urn (4.1) 
mit 


IıM=* 


? 
u—=S; Dina, 0sas ki» 
1 i=1 
Die aus sehr vielen Gliedern bestehende rechte Seite von (4.1) ist so einfach, dank den als 
Separatoren auftretenden, sinngemäße Zusammenfassungen anregenden Werten von /,. Wir 
haben etwa für s=4, r=4 


pr +) (a,0,aza,) = 2{l,p") (ar 0,0304) + ı [pP (a) pP (aa) +++] 
+1, [p (a,@,) P%(as)+- +}: 
Da die Rekursionsformeln, in denen die Werte /, eine so entscheidende Rolle spielen, biolo- 
gische Tatsachen formulieren, mag es auch aus diesem Grunde erlaubt sein rückzuschließen, 
daß unsere Segregations Wahrscheinlichkeiten nicht nur mathematisch, sondern auch biologisch 
sinnvoll sind. 

Zur Anwendung von (4.1) von Generation zu Generation brauchen wir noch Rekur- 
sionsformeln für die auftretenden Marginal W.en. Diese sind hier nicht trivial — wie 
in Abn. 2 —; trotzdem wollen wir in dieser Übersicht nicht darauf eingehen. Die Rekursions- 
formeln bilden ein System von quadratischen Differenzengleichungen, deren konstante Koeffi- 
zienten die Werte /, sind. Ein Verfahren zur Lösung, für jedes s, wurde angegeben (G. (c)), 
ebenso, für jedes s, die Konvergenzgeschwindigkeit gegen den Gleichgewichtszustand. Haldane 
hat die Rekursionsformel und explizite Lösung für s—=2 gegeben und die Konvergenzgeschwin- 
digkeit für s=3 festgestellt. 

Wenden wir uns nun der Frage stabiler Verteilungen zu. Zunächst, genau wie beim 
Diploid, bleibt die Gameten-W. jedes einzelnen Gens dauernd ungeändert 


p”)(a;) = pP” (a,), Gisele, Sr Bu re 
Ferner: Wenn komplette Koppelung der väterlichen und somit der mütterlichen Gene 
vorliegt, d.h. „=1,=}, so ist wieder 


pm(ah,... ar) = DOG. 5 Qr);, ‚0. Bart see Er ee ee 
Ist aber, im „allgemeinen Fall“ W<4, so gilt: 
lim pm(ah... ar) = (p%(a,))“r .. (PMa,))r = 2 nl), 


ImGleichgewichtszustand sind die einzelnen Allele unabhängig verteilt. 
Haldane hat, für ,=1/S, bewiesen, daß die Verteilung rechts in (4.4) stabil ist, d.h. sich 
unter dieser Segregation reproduziert; doch diskutiert er nicht — wann so eine Verteilung 
angenommen wird. — Es ist in der Tat von vorneherein plausibel, daß bei Panmixie ein Gleich- 
gewicht, wenn überhaupt vorhanden, nur in unabhängiger Anordnung der Allele bestehen 
kann. 


5. Weitere Fragestellungen für Autopolyploids. Chromosom 
Segregation ist ein Spezialfall von Chromatid Segregation, die eigentlich heute 
als zutreffend angesehen wird®). Immerhin stellt erstere, wie wir sehen werden, eine 
Annäherung an die letztere weit kompliziertere dar. Wir wollen nun diese neue Form der 
Segregation erklären. Für m=1 besteht ein Gamet aus s>1 Chromosomen. Der Zygote 
besitzt zwei Gruppen von je s Chromosomen, jedoch 4s ‚‚Chromatiden‘“, je zwei davon zu- 
sammengehörig, d.h. zum gleichen Chromosom gehörig (‚‚sister chromatids‘). Bei der Gameten- 


bildung werden s der 4s Chromatide weitergegeben, was auf S— ) Arten geschehen kann. 


Eine solche Gruppe kann von s—o verschiedenen Chromosomen stammen, wobei für s=2u, 
oder s=24+1 für o die Werte 0, 1,... win Betracht kommen; ein solcher Gamet enthält 
dann o Paare ‚‚zusammengehöriger‘“ Chromatide. Nennen wir der Kürze halber einen Gameten 


„normal“, wenn 0=0, so gibt es 8,= 2%. (7) normale, und S—S, solche, die mindestens 


ein Paar zusammengehöriger Chromatide aufweisen. Wenn und nur wenn die W. aller nicht 
normalen Gameten gleich Null ist, reduziert sich Chromatid- auf Chromosom-Segregation. 


8) Es wurde schon früher bemerkt, daß sich zeigen läßt, daß für Draade m = beliebig, =, _ beide 
auf dasselbe herauskommen, 


@ 
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R. A. Fisher und K. Mather°) definieren eine S. V., die gemäß den (u +1) möglichen Werten 
von o ebensoviele Parameter enthält. 

In unserer Definition der S.V. (s. G. (f)) tragen wir außerdem, aus denselben Gründen 
wie bisher, der Proportion mütterlicher und väterlicher Gene im Gameten Rechnung. Eine 
solche S.V. enthält die allgemeine Chromosom-Segregations-Verteilung der zwei vorigen 
Abschnitte als Spezialfall, sowie eine von Haldane (s. Fußnote 7) angenommene Gleichver- 
teilung der $ Möglichkeiten (Zufalls-Chromatid-Segregation). Wir bezeichnen also mit u 
die W. daß eine bestimmte Gruppe von & mütterlichen und s—« väterlichen 
Genen weitergegeben werde, wobei der entstehende Gamet genau o Paare 
zusammengehöriger Chromatide enthält. (o=0,1,.... u)'%). Es seien alle «2 für 
dasselbe & und o einander gleich, und außerdem u —u_,. Setzen wir uWy—=u, so hat man 
Zutstürrg— Zlundis—3: 


s=2: 4w+4u+%+2w+wW=4.2!(wtr2u)+4w=1, | 

s=3: 16(w+9u)+2@aw+3w)=1 J 

Es ergeben sich nun hier Rekursionsformeln von ganz anderer Art als vorher. Man hat 
z.B. für s=3 für einen Homozygoten: 

p"+V(a}) = 16 uop")(ar) + 144 u, play) p(a) +24 wp") (at) +36 u [pla)]’ - - (5.2) 

und ähnliche Formeln für die Typen (a7 a,) oder (a,@,@,). Ihre Struktur kann wohl analysiert 


und verstanden werden; auch Integration der Rekursionsformeln, für s=2 und s=3 wurde 
durchgeführt und schließlich das Grenzverhalten für n— o festgestellt. Z.B. die limes-Formel 


für den Fall der Gleichverteilung (alle 8 Wahrscheinlichkeiten in der S.V. einander gleich): 


e651). 


72 15 
u a Pplaypla) 777 Plalpla) 6-3) 
zeigt eine von den früher gefundenen, ganz verschiedene Struktur. Eingehendere Diskussion 
würde uns zu weit führen. 

Schließlich habe ich (G. (e)) für Autopolyploide die allgemeine Aufgabe der Ver- 
erbung von m gekoppelten Charakteren, 2s Chromosome, r Allele, unter 
Chromosom Segregation durchgeführt. Dabei müssen sich einerseits, für s=1, die 
Problemstellungen und Resultate von Abn. 2 anderseits für m=1 die von Abn. 3 und 4 
ergeben. Hier lag eine interessante Arbeit von D. De Winton und Haldane!!) vor, die für 
Tetraploids (d.h. s=2) der primula sinensis für den Fall m=2 eine durch Beobachtungen 
nahegelegte Koppelungstheorie durchführen. Unser Gedankengang ist ähnlich wie bisher: 
Zunächst ist ein Gamete natürlich durch Angabe von m loci und s aus den r möglichen Allelen 
gewählte Gene gegeben (das wären ohne alle ‚‚Symmetriebeziehungen r°” verschiedene 
Typen). Es ist nun ein Unterschied, ob man bloß die m-Gruppen von s beliebig angeordneten 
Allelen gibt (zu jedem locus s Chromosome) oder s Gruppen von je m ‚‚Vektorkomponenten“ 
(zu jedem Chromosom die m Charaktere). Die in der zweiten Auffassung zu unterschei- 
denden Gameten (a,b,, a,b,) #£ (a, b,, 4,b,) sind in der ersten, beide durch (a, a,, b, b,)gegeben. 
r+ s—1\]” m s—1 

8 E 
Beim Zygoten kommt außerdem noch die von uns, hier wie auch früher, bejahte Frage dazu, 
ob mütterliches und väterliches Erbgut aufzuzeigen ist oder nicht. Unsere S.V., die ja stets 
den Angelpunkt der Theorie bildet, ist auf der ersten Auffassung basiert (was mir heute nicht 
mehr als die unbedingt richtige Wahl erscheint)!?). Es sei 


(dl (&) ac. ge: =) (te. we ı..9..64 


die W. dafür, daß innerhalb der von einem Elter abgegebenen Gruppe von s 
zum iten Jocus gehörenden Werten sich eg; „mütterliche“und s—e; „väterliche“ 
befinden, wobei i=l,... m, 5—=0,1,...s. Es ist klar, daß sich diese Verteilung für 
s—]1, resp. für m —1 auf die früher eingeführten S. V.en reduziert. 


9) „The inheritance of style length in Lythrum Salicaria“. Ann. Eugenics 12 (1943), p. 1—23. Siehe 
auch weiter Fisher: Ann. Eugenics 11 (1941), p. 31—39 und 12 (1944), p. 169-1711. 

10) Ob die Berücksichtigung der Proportion mütterlicher und väterlicher Anteile biologisch gerecht- 
fertigt ist, kann nur durch Beobachtung entschieden werden. Wenn nicht, so kann man, nachdem die ein- 
‚geführten Parameter ihre mathematische Schuldigkeit getan haben, sie leicht beseitigen, indem man bei 
festem o alle u? für verschiedene & einander gleichsetzt, was die obige Fisher Mathersche Auffassung ergibt. 

' 11) „Linkage in the tetraploid primula sinensis.“ J. Genetics 24 (1931), p. 1—44. 
12) Sie enthält De Winton und Haldanes fürs=m=r=2 definierte 8. V. nicht als Spezialfall. 


Im ersten Fall sind R = ‚ im zweiten Fall R, -( Gametentypen. 
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Wir haben sodann die Marginalverteilungen für die Gametenverteilung, sowie für die 
S.V. zu erklären und gewinnen Rekursionsformeln für die Gametenverteilung und ihre 
Marginalverteilungen. Die Ergebnisse sind außerordentlich durchsichtig. Wir betrachten als 
einfaches Beispiel den Fall, daß nur m verschiedene Allele (je eines per locus) auftreten. Dann 
ist 


DEU ASS em) (2)- (2) mc ... Am) pm (Asa... Am) 
Ey! 


0 & Em, 

(5.5). 

Die Interpretation ist einfach: l(e,...&m) asia) ist die W., daß bez. des ten Jocus eg; 
al m 

mütterliche, s—e; väterliche Gene abgegeben werden (p”)(A&, ... A,”) ist die (marginale) W. 


eines mütterlichen Gameten mit diesen Genen, und analog für pm(A}7#, ... Am). Unter 


Panmixie gibt die Summe dieser Produkte die gewünschte W. Die Einfachheit, die sich hier, 
bei allgemeinem r, m, s noch findet, geht bei Chromatid Segregation verloren. 


Wir geben schließlich, zum Fall beliebiger Allele zurückkehrend, den ‚‚Hauptfall‘ des 
Grenzwertsatzes an. Ähnlich wie in Abn. 4 kommt es hier auf /,(0) an, d.h. auf die aus 
K(&1;...&m) abgeleitete Marginalverteilung erster Ordnung. Wenn 1,(0)=4, oder 
;—=1—21;(0) = 0, so besteht bez. des «" Jocus complette Koppelung (und die betref- 
fende Marginalverteilung bleibt unverändert erhalten); andernfalls gilt: Für allgemeines m, 
r, skonvergiert die Verteilung der Gameten gegen das Produkt der ms Aus- 
gangsverteilungen der einzelnen Gene, wenn und nur wenn die m W.en c; posi- 
tivsind. Das Resultat ist von allen in die Segregations-Verteilung eingehenden Parametern 
unabhängig. 

Wir haben hier eine Gruppe von Problemen aus dem Gebiet der Populationsgenetik 
behandelt, die sich auf Basis der Mendelschen Theorie, unter Benutzung eines rationellen 
Wahrscheinlichkeitsbegriffes, mathematisch beschreiben lassen. Natürlich sind die ein- 
geführten Begriffe und Hypothesen nur als mathematische Modelle der Wirklichkeit zu 
betrachten, und zwar als ,‚Rahmenmodelle‘, derart daß man durch geeignete Verfügung über 
die vorhandenen Parameter mehr konkrete, spezifischen Situationen angepaßte Darstellungen 
erhält. Anderseits mag es erwähnenswert scheinen, daß manche Resultate unabhängig von 
allen Parameterwerten bestehen. — Zwischen der Mathematik und einem stets wachsenden 
Kreis von Wissenschaften besteht das Verhältnis gegenseitiger Anregung, Hilfe und Berei- 
cherung. Es ist hier auf einem beschränkten (doch erweiterungsfähigen) Gebiet der Versuch 
gemacht, die mathematischen Werkzeuge der Abstraktion und Verallgemeinerung in den Dienst 
unseres Strebens nach immer besserem Verstehen zu stellen. 
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Über die Theorie der dünnen, schwach gebogenen Platten 


Herrn R. von Mises zum 70. Geburtstag 
Von @eorg Hamel in Landshut 


Die verdienstvolle Arbeit von Herrn Manfred Schäfer in dieser Zeitschrift (Bd. 32, 
Heft 6), die uns mit den grundlegenden Arbeiten von Herrn Erich Reissner über eine 
verfeinerte Theorie bekannt macht, veranlaßt mich, einiges Grundsätzliche zur Theorie der 
dünnen Platten zu sagen. 

Beide Theorien, die klassische von Kirchhoff sowohl als auch die verfeinerte von 
Reissner enthalten das Unbefriedigende, daß außer genauen Voraussetzungen auch Vernach- 
lässigungen gemacht werden, bedingt durch den Übergang vom Raum zur Fläche, deren 
Tragweite und Berechtigung schwer zu übersehen ist, z.B. das Fortlassen einer ganzen 


OÖ 
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Gleichung (5c bei Herrn Schäfer). In der klassischen Theorie ist die unangenehme Folge 
wohl bekannt: man kann nur zwei unabhängige Randbedingungen befriedigen, während drei 
bestehen. 

Dieser Übelstand wird von der verfeinerten Theorie überwunden. Aber der erreichte Grad 
der Annäherung erscheint noch prekärer. Man beachte nur, daß in der Gleichung (II), aber 
auch schon in den Gleichungen (13d) u. (13e) des Herrn Schäfer ein kleiner Parameter bei 
den höchsten Ableitungen steht, was bekanntlich für eine allgemeine Theorie höchst mißlich ist. 

Demgegenüber möchte ich darauf aufmerksam machen, daß es eine in sich geschlossene 
Theorie der dünnen Platten gibt, die rein zweidimensional arbeitet und auf Grund präziser 
Voraussetzungen völlig streng vorgeht. Sie ähnelt zwar in den Ergebnissen der klassischen 
Theorie, vermeidet aber die genannte Schwierigkeit; sie kann mühelos ohne ‚künstliche 
Überbrückung“ drei Randbedingungen befriedigen. 

Zu dieser Schwierigkeit möchte ich zunächst noch etwas sagen. Die Bedingungen (12d) 
und (12e) der klassischen Theorie, nämlich 

ow ow 
Ge di Zn 
dürfen nicht als Folgerungen der dynamischen Gleichungen (5) angesehen werden, sondern 
als Folgerungen der Voraussetzung, daß die Normalen Normalen bleiben. Sie sind also kine- 
matischer Natur. Es wäre demnach unbillig, für sie dynamische Randbedingungen wie (17a) 
zu fordern. Die Größen ®, und ®, kommen auch in der strengen Theorie gar nicht vor. 

Diese strenge, rein zweidimensionale Theorie habe ich in Kap. III, $4 meiner Theore- 
tischen Mechanik dargestellt. Eine erste Form, die in den Abh. der Berl. Akad. d. Wiss. 
veröffentlicht wurde, leidet an einem Vorurteil, das mit dem Kernpunkt zusammenhängt, 
den ich gleich erwähnen werde. Ich bitte, diese Abh. nicht zu Rate ziehen zu wollen. 

Da der $4 des Buches sehr umfangreich, sehr formelreich und doch noch sehr knapp 
ist, gestatte man mir, ihn auszugsweise hier hinzuschreiben, wobei ich einige Ergänzungen 
hinzufügen werde, die ich im Buche dem Fleiße des Lesers überlassen habe. Ich werde aller- 
dings die Beweise fortlassen, diese muß man im Buche nachlesen. 


Zusammenstellung der Ergebnisse: 


In der klassischen Theorie gelten in der exakten dagegen 
die Differentialgleichungen 
+ tp=0  (a)| Fr fr 4 p=0 
n 0 (4b) er Q=0 
+ a rt. 


Hier sehen wir gleich den Kern des Unterschiedes: es gibt keinen Grund 7,+ 7,=0 anzu- 
nehmen. (In meiner ersten Abh. war ich hier selber noch befangen.) T,+ T,#0 zu beachten, 
ist keine willkürliche Annahme, sondern exakte Folgerung aus der klassischen Lagrange- 
schen Mechanik. T, und T, sind die Torsionsmomente in den Schnitten 2=const bzw. 
y= const. 

In der exakten Theorie gibt es wohl eine wesensverwandte Beziehung statt der nicht- 
geltenden. Für beliebige Schalen heißt sie 
D’ 
A 
(&, ß beliebige Koordinaten auf der Fläche, E,F,G und D, D’, D’ die sechs Fundamental- 
größen, 


F 


(«=D F—D9); 6, =, (DF—D’B) M—-EG-R, 


S2, Sz die tangentialen Schubkräfte in den Schnitten, Z,, Z, die Zugkräfte.) 
In unserem Sonderfall wird daraus: 
ow 02w 


0° 
+ eo), IT. mudy Oy? =0, 


also in erster Näherung S, + 8, =0, was uns nicht weiter berührt. 
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Wenn wir die Kirchhoffschen Annahmen zu Grunde legen (s. Schäfer, S.163 unten) 
und noch von der Querkontraktion absehen, was eine unwesentliche Vereinfachung bedeutet, 
bekommen wir die 

drei Randbedingungen: 


2 | 2 2 y\2 
WA; hs ie eich 5 Adach 2Blo + —e (er) 


0x 0y ds ds 022 ° oya/l\ds 
: w rw ee] 22)‘ 
ll N —— he ee _ -_ D—- Ma 
z ;»[42, Her Her 0y%/I\ds (A) 


A, B, sind Materialkonstante, M das Besungmonen am Rande, 


16 oAwdy | In karl Ra da dy 
jet 1° U 8 LER 
D el an da Fon. de EN " zaral au Dya) ds ds 


$ 
9?Ww (2) je a jr 
09x 0y \\ds [a5 

Q die Normalspannung am Rande, 7 das Torsionsmoment. 
Diese Gleichung fehlt im Buche, sie wurde als Aufgabe gestellt. Gleichung (A) findet 


sich S.166 oben. Integriert man nach s, so ist die Eindeutigkeit von selbst gewahrt, wie man 
leicht nachrechnet. Endlich noch 


ya ET, et: en ESBAIEH | 
“ "(Er a) Gr 9 3} i Al 002 Hoy/dsds 0x09y\\ds ds | (er 


Wie in der klassischen Theorie besteht auch hier eine Gleichung der Form 


1 


- (B) 


const Adw=p9, 


- Ka). 4 : ; Ä > Lea: 
wo die Konstante 3 h®B ist. (Bei Bestimmung der Konstanten spielt die Vernachlässigung 


der Querkontraktion eine Rolle; bei einer rein zweidimensionalen Theorie wird man 4, B, 
als Materialkonstante ansehen, deren Wahl zur Anpassung an die Erfahrung dient.) 
Ist diese Differentialgleichung mit Randbedingungen für w integriert, so kann man 


I IT 
aus (A), (B), (C) die Größen M, u und — Fr berechnen. Man kann aber auch M, Q und T, 
nur den Bedingungen der allgemeinen Mechanik entsprechend, vorschreiben, bekommt zwei 
Randbedingungen für w gemäß (A) und (B) und kann dann noch 7,--T,, berechnen, das im 
allgemeinen nicht null sein wird. 


Gedankengang des Beweises 
Die Fläche sei durch den Vektor 
r=1r(0,ß) 
gegeben, ein Nachbarpunkt auf der Normalen » durch 
n=r+ny yIsh. 
Die räumlich verteilten, eingeprägten Kräfte seien 


xdV =xAdadßdy 
so daß 
h 


—h 
wird. Ihre virtuelle Arbeit ist 
64, = [[[#(5°+ y On) Adadpdy—2h || 67 Adadß 
h 
da [ydy = 0 ist. 
en 


Nach den Prinzipien der Lagrangeschen Mechanik, insbesondere dem Befreiungsprinzip 
hat man zunächst alles so starr wie möglich zu nehmen, d.h. 


3E=0 (1), SF=0 (2), 8G=0 (3) (Unausdehnbarkeit) 
oeD=0 (6), $8D’=0 (7), 58D”=0 (8) (Unverbiegbarkeit). 
Weiterhin nehmen wir mit Bernoulli an, daß die Normalen solche bleiben, d.h.: 
EMO. 5... (4), undarıdln 20 Bee 


oO 
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wo € = FR Fr 58 sind, und weiter wollen wir noch, was unwesentlich ist, von Querkon- 
traktion absehen, d.h.: 
Be) = DEE URL... Bern (9) 
annehmen. Schließlich soll der Rand ganz festgehalten werden, d.h. an ihm 
GE er l103] undald 0 .. .... . ..% (106) 


gelten, wo 68 eine ®irtuelle Verdrehung ist. 
Man hat nun diese Bedingungen (1) bis (10b) mit Lagrangeschen Parametern 


1 1 


hu» Aa, IE Ads, Mds 


ar nn 
2;: Ass As; A h—=A + 8,, de; 5) 
4 


zu multiplizieren, zu integrieren und zu ö4, hinzuzufügen. A, muß deshalb in zwei Parameter 
zerspalten werden, weil 6D’ zwei formal verschiedene Gestalten haben kann, nämlich einmal 


eönor cm oör 


Ar 08 2a oß 
aber auch 
oönder Omdäör 
D—= — E 
) 5 oB da OB da ' 
A, gehört zur ersten, A5 zur zweiten. (Das hatte ich bei meiner ersten Publikation übersehen.) 
Setzt man 


u —- h =2u, 


so gehört — u als Lagrangescher Parameter zu der rein mathematischen Relation 


(7) — JEDE ER. BR 


Man hat nun den um die Glieder mit den Parametern vermehrten Ausdruck der virtuellen 
Arbeit durch partielle Integration so umzuformen, daß nur noch Glieder mit ör, ön, 6# vor- 
kommen und diese einzeln null zu setzen und zwar sowohl bei den Flächenintegralen wie bei 
den Randintegralen. Man erhält so zwei vektorielle Feldgleichungen und zwei vektorielle 
Randbedingungen, die ich herschreibe: 


Sure merke 
RN BL AND-aD-—AR)— in) 
* Tess 5a 2 OB 7, A 8 5% 
1) on 0 0 (., OR () on 
Y +8; ET en (a tl) = re) 
MEHR HRD + BE oz UI nleh Bel LI), 
Randbedingungen: 
TERN, on „„onld RN... d 
SE Re Rh ee -45]% —+4=0 (II) 
wi d : 
nxı— (Ag € LITER IR: en near. (IV). 
Dynamische Bedeutung der Parameter: 
Mit 
— or or . o > 
M'== MZ + Tz (v äußere Normale in der Fläche) 
sind 
@ , E E F R F 
= M. A a a: h=T, 74 "=—-R—-— Ms 


el: 
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Mit 
A =! RN 
0v os 
,= —Vad,:; = — YEQ,. 


Für Z,, Za, 8%, Ss ergeben sich vier Gleichungen, aus denen sich A,, A,, A; ausrechnen lassen, 
was uns hier nicht besonders interessiert, A, ergibt sich doppelt, was die oben genannte Gleichung 
zwischen 8%, Sg, Zg, Ze, Ma, Ma; Ta, Tg zur Folge hat. 

Nun hat man zu befreien, d.h. die Bindungen (1), (2), (3), (6% (7), (8) aufzuheben, 
die Fläche also dehnbar und biegsam zu machen. Damit werden A,, Ag, Ag; Ag, Ar, Ag eingeprägte 
Kräfte, nicht dagegen A,, A,, wenn wir an der Bernoullischen Annahme festhalten, daß die 
Normalen Normalen bleiben und auch nicht A,, wenn es bei Vernachlässigung der Querkon- 
traktion bleiben soll. Besonders zu beachten ist, daß zwar A, — A, +47 eingeprägt wird, aber 


1 
nicht u => (45 — 47), das zu der erwähnten rein mathematischen Beziehung (11) gehört. 
Bei Annahme vollkommener Elastizität wird demnach 


1 


De N LIT 55, =79Dr KDD 


zu setzen sein, wo öl fydadß die virtuelle innere Arbeit ist. Setzen wir noch Isotropie 
voraus, so kann y nur von den beiden Invarian*ten abhängen, die wir so erhalten: 


Aus Ab: ÖO 


N=r+Nny 
(dr,)®—= E*da® +2 F*dadp + G* dB? + dy? 


berechnet man 
mit 2 
2 
on on 
*_ ee 4 2 
a 24 (r® Ar =), +)» 
G=G-+2 4 (2 y: 
De Truck 
und daraus die absoluten Invarianten: 
E* @* — F*? d E*05 + EIG*— 2F*F% 
Bam BG re 
wo die mit dem Index O versehenen Größen die im Urzustand sind. Mit 


Br—Bt+er, PR 4fr, GG 4g* 


sind auch 
e* g* — f*? e* + g* ES 
, ES 2 ET Bar Ter 
IS po m Wi B=-—a—re 


Invarianten und daher bei weiterer Annahme eines linearen Gesetzes nach Hooke: 


h h 
v=Jpdy= (AI + BI2)dy 
A, B, Materialkonstante. 

Es ist nun leicht, im Sonderfall der Platte mit der Benutzung von a—=z, ß=y und der 
Kirchhoffschen Annahme, daß keine Verschiebungen u, v in der Platte, sondern nur eine 
solche w senkrecht zur Grundebene stattfindet, die weitere Rechnung unter den üblichen 
Vernachlässigungen bei kleinem w durchzuführen. Es ergeben sich 


h=0, ,=0, ,=0 bei Ad —=0, Aj7—=0 am Rande, 
o°w =) 
EIERN N \ a 
47 uee u ei 
8 0 
her ren 


DW 08 
ee Sujarı + lt]. 
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und daraus alles weitere. Aus den zwei vektoriellen Feldgleichungen sind die Reaktionen 
Ay, As, u zu eliminieren, so daß, weil zwei Gleichungen wegen A,—=0, A,=0, /;,—0 von selbst 
erfüllt sind, nur noch eine für w übrig bleibt, nämlich das oben angebene 


= BAdu— ». 


Von den Randbedingungen bleiben noch drei stehen, die in zwei für w und eine für das noch 
unbekannte u umgeformt werden können, wie oben angegeben. 


Elementarer Existenzbeweis für die Strömung in der laminaren 
Grenzschicht zur Potentialströmung U=u,2” mit m>0 bei 
Absaugen und Ausblasen. 


Herrn Professor Dr. Richard von Mises zum 0. Geburtstag verehrungsvoll zugeeignet 
Von Rudolf Iglisch in Braunschweig 
1. Einleitung: Die Formel 
Be EEE 07 Eee) A. ru. (1) 


und konstantem u, liefert (vgl. z.B. [1]) die &-Komponente der Geschwindigkeit bei einer 
Potentialströmung symmetrisch um einen ebenen unendlich ausgedehnten Keil vom Öffnungs- 
winkel Br mit 


2m ß 
Br Se u! . . . . . . . . . . . . 2 
wenn die z-Achse in eine Keilkante mit der Spitze als Ursprung gelegt wird, die y-Achse 


senkrecht dazu in die Strömung hinein. Nach V. M. Falkner und S. W. Skan [2] 
werden die Prandtlschen Grenzschichtgleichungen erfüllt durch den folgenden Ansatz für die 


Stromfunktion: 
\ 2 m+17/uam-1 
— | ————E m gr PER SE 


(v = kinematische Zähigkeit), wenn f(n) der gewöhnlichen Differentialgleichung! 


a Fantasien) 
genügt. Für die Geschwindigkeitskomponenten findet man: 
| Be ee 1 
u=marf(n): = Yuar ( 1:7 u) BER (0): 
Zu der Differentialgleichung (4) nehmen wir noch die Randbedingungen 
ee et a, ei erellejene « (6) 


hinzu, von denen die beiden ersten vermöge (5) aussagen, daß die Flüssigkeit keine Geschwin- 
digkeitskomponente längs der Wand besitzt, wohl aber eine ganz bestimmte von & abhängige 
Geschwindigkeitsverteilung senkrecht zur Wand, und zwar bedeutet Ü >O Absaugen, Ü<O 
Ausblasen; bei anderen Absauge- oder Ausblasegesetzen läßt sich das Problem im allgemeinen 
nicht auf ein eindimensionales reduzieren. Die dritte Bedingung (6) sichert wie üblich den 
Anschluß der Grenzschichtströmung an die Potentialströmung. Man erhält Strömungsvor- 
gänge ohne Rückströmung, wenn man noch die Zusatzvoraussetzung 


an.) nad tl?) 
hinzunimmt. 

Das Problem (4), (6) enthält die beiden reellen Parameter ß und (, von denen im folgenden 
ß=0!) angenommen werden soll. Für die Wertepaare ß=0, C=0 und f=1, C=0 hat 
H. Weyl[3, I] 1941 die Existenz einer Lösung, die der Zusatzbedingung (7) genügt, mit Hilfe 
eines Verfahrens der sukzessiven Approximation gezeigt. Im Fall 20, Ü=0 wurde die 
Existenz einer analogen Lösung 1942 in [3, II] bewiesen, und zwar unter Heranziehung von 


1) Auf die Einschränkung ß < 2 kann verzichtet werden. 
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Fixpunktsätzen über Funktionaltransformationen, wie sie zuerst von G. D. Birkhoff und 
O0. P. Kellogg [4] angegeben und von J. Schauder und J. Leray [5] in großer Abstraktheit 
formuliert wurden. Eine ausgedehntere Darstellung fanden die Beweismethoden 1942 in [6]. 
Außer in den Fällen ß=0 und ß=1 wird die Frage der Eindeutigkeit der Lösung nicht 
behandelt. 

Ohne infolge der Zeitumstände die Arbeiten von H. Weyl kennen zu können, bewies 
ich 1943 Existenz und Eindeutigkeit der Lösung von (4), (6) unter der Zusatzbedingung (7) 
für die Wertepaare ß—=1, (beliebig [7]. Daß ohne die Zusatzbedingung (7) die Eindeutigkeit 
(wenigstens für genügend große positive C') nicht mehr gilt, wurde in [8] gezeigt. Beide Arbeiten 
benutzen nur elementare gestaltliche Aussagen über die Lösungen. Mit anderen Hilfsmitteln 
wurde von R. Iglisch und D. Grohne [9] gezeigt, daß im Falle $=0 das Problem (4), (6) 
unter der Nebenbedingung (7) nur für © 2—1,2385... lösbar ist. Ziel der vorliegenden Note 
ist es, die in [7] und [8] entwickelten Schlußweisen noch durchsichtiger zu gestalten und auf 
das Problem (4), (6), (7) für # z0 und alle Werte von ( anzuwenden. Da der Eindeutigkeits- 
beweis — vor allem für die Wertepaare (ß, ©) mit ß<4 und O<0 — langwierigere Über- 
legungen erfordert, muß er aus Raummangel einer späteren Veröffentlichung vorbehalten 
bleiben. Diese Note enthält den elementaren Existenzbeweis für > 0 und alle Werte von (; 
da nach [9] im Fall #=0 nicht für alle Werte von (' Lösungen existieren, wird man nicht 
erwarten können, daß die elementaren Methoden ohne spezielle Anpassung diesen Fall mit- 
liefern. 

Es sei gleich hier erwähnt, daß das Problem (4), (6), (7) mit 


FN)=f(n; f’=FF, f"= (Fi) F— FrLIr me ea 
(- bedeutet Ableitung nach f) gleichbedeutend ist mit der Lösung der Ditferentialgleichung 
F(FB+/FE+BA—P)—=0. 2 
unter den Bedingungen 
FO=0,=-H8)= 1; Hifye> 0 für f > Ges De 
Aus F(f) erhält man mittels 
1 [6 an) FT 
den Zusammenhang zwischen f und n. 
2. Hilfssatz 1: Die Lösung F(f) von (9) unter der Anfangsbedingung 
FOSU.LFO FO Seel 


Akira bei >20 und genügend großem a für gewisse Werte f>( Funktionswerte an, die 
größer sind als 1. 


Beweis: Der Satz ist gleichbedeutend mit folgendem: Die Lösung von (4) unter den 


er 
(0) = 0,7 PO) = 02T 101 0 ee 
besitzt bei $>0 und genügend großem a Ableitungswerte f(n)>1?). 
Zum Beweis des letzteren Satzes unterscheiden wir zwei Fälle: 
u Ba a 


Wir betrachten nur das Intervall O<n<1. Wäre hier im Gegensatz zur ae 
stets die Annahme 


fmsı für nina re 2 
richtig, so wäre dort 
IM) zsCOFr<ZC FT. rm Er 

Nun ist nach (4): 
f’=--/f’—-Bi1—f?) s0, solange ’(n)z0istt ...: . (16), 


wie nach (14) folgt. f’’(n) nimmt also solange von f’(0)—=aanab. Daher folgt aus (15) und (16), 
solange f’(n) Z0 ist, 
KaaU Eu 162 297 2 
2) Daß die hier und im folgenden benutzten Lösungen der De Ran existieren, soweit 


sie benötigt werden, ist leicht zu sehen. Es folgt z. B. auch aus einem Satze von F. Rellich [10] (vgl. dort 
insbesondere $. 171 unten). 


OÖ 
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daraus durch Integration 
Bine |(O Ella Bine ne ano sale. (17). 
Sicher ist mithin f”’(n) 20 bis zum Wert 
a 
na 
hin. Wir betrachten fortan nur das Intervall O<nsnm si. 
Hier folgt aus (17): 


(18) 


Sn) zan—3lC+Na+Pß]r 
und somit nach (18) 
, a? 1 a? 1 a? 
2 Cm BeccHmars 
bei genügend großem a entgegen der Annahme (14). 
Eall- 11: 0.0: 
Man schreibe (16) so um: 


"= 4-9 BA) 0”. 


Nun läuft der Schluß analog wie früher. Man erhält, wenn man jetzt entsprechende Formel- 
nummern mit einem Stern versieht: 


a One Tee Re rn ne. (15*). 
Solange f’ >0 bleibt, folgt aus (19): 
lin 208,, f (0) a, 
damit 
n 
Po) = er=9 [a+ [(--pos)e-=0n di], 
N) 
also fürOsn<sI 
m) saet. 
Somit folgt aus (19) und (15*): 


ee ee ma a [ae EBln re 2 . (17%). 
Fm)=0 bis 
gt an 0 N 18r 
hin. Schließlich wie vorhin 
1 a? 
en) az! 


für genügend großes a, entgegen der Annahme (14). 

3. Hilfssatz 2: Die Lösung des Anfangswertproblems (9), (12) bleibt bei genügend 
kleinem a stets kleiner als 1, falls $>0 ist; sie steigt monoton bis zu einem Maximum <1 
an und fällt danach stets ab, solange (7) erfüllt ist. 

Beweis: Fall I: C >20. 

(9) geht mit der Abkürzung 


Be BIBBen ,.,..... (19) 
infolge 
dU dU, 
df dF 
über in 
u = BR) —p(M) 2 NE NA 1) 


wobei p(F) = fU infolge f 2 C eine positive Funktion ist, solange U >Ogilt. Integration von (20) 
liefert unter Beachtung von (12): 
De 5g® F®3 a 
ee te u ie (21). 
en) (21) 
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Da bei genügend kleinem $a®? der Ausdruck rechter Hand zwischen F=0 und F=1, wo er 
verschiedenes Vorzeichen hat, verschwindet, wird daher an einer Stelle #, zwischen O und 1 


die Funktion U= 0, mithin dort auch #=0. — Sei etwa F(f})= 0. An jeder solchen Stelle fı 
mit F=0, O<F<1 ist aber nach (9) und (8) F<0, es liegt also ein echtes Maximum vor; 
hinter fi muß daher F(f) dauernd abnehmen, solange F >0 ist. 

Falltze0z2y: 

Man schreibe (9) so um 


FFF+(+lC)) FF— OFF +BA—P)=0. 
Dann erhält man 


u = eo rest) nee ....(20*), 


solange U >0 ist, und es kommt wieder nur darauf an zu zeigen, daß U(F,)=0 sein muß 
für einen Wert F,<1. Die rechte Seite von (20*) ist für #=0 bei genügend kleinem a negativ; 
daher nimmt für kleine F-Werte U ab bis zu einer Stelle F,, an der die rechte Seite von (20*) 
verschwindet. Für O<F<SF, gilt sicher 

ICIU < IC U(0) = |Ola 
Solange kann man statt (20*) schreiben 


u - Cla—p(1—F)—p(F)—p(F) ...... . (20*#) © 


und findet zunächst 
r2/1-05=R<1. 


Durch Integration von (20**) erhält man 
LazeRas F3 
_—— — _— -— ng . . . . . . . . . . . . ” . 
5 <5-6-I0aF+B7 21”) 
Da bei genügend kleinem a die rechte Seite bei #=0 und F=F, verschiedenes Zeichen hat, 
folgt wieder die Existenz einer Stelle F, der verlangten Art. 

4. Hilfssatz 3: Es gibt im Falle $>0 keine Lösung F(f) von (9), die sich für große 
Werte von f monoton wachsend einer horizontalen Asymptote #=b mit 0<b<{1 annähert. 

Beweis: Indirekt. Für f2/%,>0 sei schon 1>b>F(f)>b—e>0. Aus (9) und (8) 
folgt, daß F>0 für f>0 und O<F<1 nur gelten kann, wenn E<O ist; also ist in unserem 
Falle sicher F<0 für f>f,. Aus (9) ergibt sich weiter für ff 

Se BR; >ß(1—b), 
PET 
pn ze Te) Sue @3 
Integration liefert 7 
—Enz F(f) >e(f— Po): 
FA)>e(f—R): 
was bei großem f zum Widerspruch führt. 

5. Der Existenzsatz: Das Problem (9), (10) besitzt im Fall $>0O mindestens eine 
Lösung F(f). — Daraus folgt: Das Problem (4), (6), (7) besitzt bei ß>0 mindestens eine 
Lösung /(n). 

Beweis: Unter den Anfangsbedingungen (12) hat (9) nach Hilfssatz 2 bei genügend 
kleinem a Lösungen, die bis zu einem Maximum <1 ansteigen und dann abfallen. Sei a* 
die obere Grenze aller a-Werte, so daß für «<a* das Maximum <1 bleibt. Ein solches a* 
muß in Hinblick auf Hilfssatz 1 existieren. Zu a* gehöre die Lösung F*(f) von (9). 

I. F*(f)kann F =1 nicht im Endlichen berühren, da F*(f)=1 die einzige Fösune, von (9) 
ist, die an einer Stelle f, bei horizontaler Tangente den Wert 1 annimmt. 

II. F*(f) kann kein Maximum besitzen für einen enalichen Wert fı. 

a) F*(f)kann keine Werte >1 annehmen wegen der stetigen u ee der Lösungen 
von (9) von den Anfangsbedingungen. 
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b) F*(f) =1 war unter I ausgeschlossen. 

c) Wäre F*(f)=d<1, so gäbe es ein f,>f,, mit F*(£)<d, etwa =d—2e. Man wähle 
a=a*+ö. Wegen der stetigen Abhängigkeit der Lösung von den Anfangsbedingungen müßte 
für das zugehörige F,(f) bei genügend kleinem ö gelten: 


d+te>F,(f)>d—e F(fp)<F*rf)te<d—e, 
was für alle #, mit genügend kleinem ö bei Beachtung von F,(C)=0 das Vorliegen eines 
Maximums <1 bedeutet gegen die Konstruktion von a*. 

III. F*(f) muß sich also monoton wachsend für f— oo einer Asymptote <1 nähern, 
die nach Hilfssatz 3 nicht <1 sein kann. Also ist F*(f) eine Lösung des Randwertproblems 
(9), (10). ; 
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Der Aufschlag von Rotationskörpern auf eine 
Wasseroberfläche 


Richard v. Mises zum 70. Geburtstag gewidmet 
Von 0. Schmieden in Darmstadt 


Der Aufschlag von Rotationskörpern auf eine freie Wasseroberfläche wird unter Voraussetzung schwacher 
Konturkrümmung untersucht und gezeigt, daß sich die Integralgleichung des Problems auf eine Abelsche 
Integralgleichung reduzieren läßt. Für einige Spezialfälle wird die Lösung in geschlossener Form angegeben 
und die auftretende Stoßkraft berechnet. 

The striking down of rotary bodies on a free water-surface under ihe supposition of a weak curvation 
of the contours is examined, and it is shown that the integral equation of the problem can be reduced to an 
Abel integral equation. For some special cases the solution im closed form is indicated and the impulsive 
force acting thereby is calculated. 

Le ricochet de corps rotatoires sur une surface libre d’eau & la supposition de courbement faible des 
contours esi examine, el il est montr& que I’Equation integrale du probl&me peut ätre reduite & une Equation 
integrale d’Abel. Pour quelques cas sp6ciau:x la solution en forme serr6e est indiqu£e, et la force de propulsion 
qui se produit est calcul£e. 

Yıaap Te Bpaımenns, TaNamımmXx Ha CBOÖOAHYIO BOAHYIO IOBePXHOCTB, HCCHENYETcH B IIPe]- 
HOJOXeHHH CıIa00U KPHBU3SHEI KOHTYP&; IOKA3bIBAeTCH, YTO UHTETPANBHOS YpaBHeHHe 3anayn 
IIPHBOAUMO K A0eJleEoMmy MHTerpalbHoMy IIpaBHeHmm. ‚(1 HEeKOTOPBIX YAcTHEIX CILyyaeB 
pemeHue NaHO B 3aMKHYTOM BU U BEIYHCHEHA BOSHUKAMINAA CHA YAapa. 


Läßt man einen Körper auf die freie Oberfläche einer Flüssigkeit aufschlagen, so ergeben 
sich bekanntlich sehr große Stoßkräfte. Um diese Kräfte wenigstens in erster Näherung zu 
erfassen, setzt man nach dem Vorgang von Wagner!) den Körper als so schwach gekrümmt 
voraus, daß die Spritzerbildung keine wesentliche Rolle spielt. Die von Wagner für den 
zweidimensionalen Fall entwickelten Ansätze lassen sich sinngemäß auf gewisse einfache drei- 
dimensionale Fälle übertragen ?), in denen auch der rotationssymmetrische Fall enthalten ist; 


ı) H. Wagner, Z. angew. Math. Mech. 12 (1932), S. 193. j R 
2) C. Schmieden, Ing.-Archiv 10 (1939), S. 1. Im folgenden mit „I“ zitiert. 
11* 
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alle für die Berechnung der Stoßkräfte charakteristischen Größen wurden aber damals nur 
in Form von Reihenentwicklungen gewonnen. Hier soll nun gezeigt werden, daß für den Fall 
der Rotationssymmetrie die beherrschende Integralgleichung vom Abelschen Typus ist und 
über eine gewöhnliche Dgl. 1. Ordnung in einigen Fällen geschlossen gelöst werden kann. 


Anstatt den Rotationskörper auf die Wasseroberfläche aufschlagen zu lassen, kann man 
sich auch die ursprünglich ebene Flüssigkeitsoberfläche gegen den ruhenden Körper bewegt 
denken, wobei die Geschwindigkeit V(t) im Unendlichen als Funktion der Zeit t gegeben sei. 
Bei diesem Stoßvorgang ist die Geschwindigkeit an der Oberfläche wesentlich vertikal ge- 
richtet, während sie am Körper selbst tangential zu ihm verläuft. Ist der Körper nur schwach 
gekrümmt (durchgehend kleine Tangentenneigung zur Oberfläche), so kann diese Strömung 
in jedem Augenblick angenähert ersetzt werden durch die Strömung um eine Kreisplatte, 
die genau die Größe der momentanen Stoßfläche besitzt. Die Geschwindigkeiten an der freien 
Oberfläche sind dann mit demselben Grade der Annäherung zu ersetzen durch die Vertikal- 
geschwindigkeiten in der Symmetrieebene der Plattenströmung. 


Die Strömung um eine in Richtung ihrer Normalen angeströmte Kreisplatte ist bekannt; 


ist cihr Radius, so gilt für einen Punkt R>cin der Symmetrieebene, daß die Normalgeschwin- 
digkeit v„(t) wird (I, 48): 


12 [6 SE 
=] 14 A are) | n. KNE Ta REREE  CE 


Das Wasserteilchen an der Stelle R erreicht zur Zeit t eine Erhebung über dem Nullniveau 
n(R,t): 


c 


(RW = [mo a a 


Es erreicht offenbar die Oberfläche des Körpers, wenn die momentane Stoßfläche die Grundriß- 
projektion des Teilchens erreicht, also für c= R. Setzt man noch 


vu) ee 


wobei « als Quotient zweier Geschwindigkeiten dimensionslos ist, so erhält man für «(c) 
folgende Integralgleichung: 
R 


n(R)= j «ol1+z Fe 


0 


C 


en) | 


Hierin ist n(R) die Gleichung des Meridianschnittes des Rotationskörpers. Gl]. (4) läßt sich 
nun in eine Abelsche Integralgleichung verwandeln; dazu wird eine neue Funktion y(c) 
eingeführt, die definiert ist durch 


[wo de = ey; RE Rn, 


Beachtet man, daß y(0) endlich sein muß, so geht (4) durch partielle Integration über in: 


d 
n(R) = 2 [roten et 6) 
Wird hier RP= x; ® = gesetzt, so folgt mit 
M=oa 90-0 ed=2e, 
2 f dy] d 
vo ere 0 


und mit der Abkürzung 


EArIOEW HJ 17 ee) 


Ö 


© 


De 
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die Abelsche Integralgleichung: 


a. 
Ihre Lösung lautet bekanntlich 


n:f(&) = IE: Ber EINER (10), 


Setzt man statt x und & wieder R und c ein, so folgt 


> dn 
e Zul 
2y()+e7 RT Bere. (11), 


R: 


Das ist eine lineare Dgl. 1. O., deren Integrationskonstante verschwinden muß, da die Lösung 
der homogenen Dgl. proportional c”? ist. Für einige spezielle Fälle soll jetzt der Zusammen- 
hang zwischen n(R) und «(c) bestimmt werden. 


TsRall: 
n(R)=k,2;. v21. 


Damit wird (11): 


e I 
v—l » Kny ee! ne | 
2 yl)+ 0? — vhy S SAETSUE | Ei Kl p)TUR ap. 


Ve IR? 2 
) 


1 el 
Be Re (3) 76) 
v/2— en. 2— u 

Te 2, DT 

0 De) 


Da aber 


ist, so wird 


2A ylc)+ er 10m; y= 


2 r ver N 
Hieraus folgt dann 
Ay 1 vA, v—1 ; 
= ea  — en 12 r 
vo=T7et Warte =.Te (12) 


In (12) steckt speziell der Fall des geraden Kreiskegels; ist & der halbe Öffnungswinkel, so gilt 


MEILE, u()= Teig. 


2. Fall: Die Kugel; wird der Kugelradius gleich eins gesetzt, so gilt 
— — - 1 
„(RB=1-fi—-R; oa)=1-Y1- P@)= BE ; 
2yı —% 
also nach (10): 
& e 
1 da 1, 1+YE& dv 1, 1+e 
= — | ee = — In ——; 2Y()+c—-=—1in ———. 
af) er er) Anl He) de 2 1-ec 
0 


Durch Integration dieser Dgl. erhält man 


1. 1-@ 1+e, rer imo 1: . 
Mh a de "rat Meer in, Far? . (13). 
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3. Fall: Die Spitze des Körpers ist nach einer Sinuslinie geformt: 


n—= 4Aarcsin R, (ERST 
Nach (11) wird hier 


e 1 
2 em A 2 
de Y( — R)(@— Ro) Yd-a)(1- ea) 
0 0 
Durch Integration erhält man 
E(c) — (1— e)K(ec 
Eu 2) Re 


Darin sind E(c) und K(c) vollständig elliptische Normalintegrale. 


Berechnung der Stoßkräfte. Das Potential der senkrecht zu ihrer Ebene ange- 
strömten Kreisplatte lautet (I. (1) für5b=e) 


2 (ec e\| 
B=Vall+ (7 arte) Pe 
Barlır ya ei 
Hierin ist & die Normalenrichtung und c der Radius der Platte, während A mit den Koor- 
dinaten x, R in einem Meridianschnitt verknüpft ist durch 


BOT re A 


K(c) — 
u(c) = A ee 


x m 
7 a ee 
In der Plattenebene ist außerhalb der Kreisplatte das Potential Null; auf der Platte selbst Ö 
sind x und A Null, aber derart, daß 
f x 1 R2 
im— =1— —, 
x, 0 4 e° 
Damit wird auf der Platte das Potential — und zwar auch der Absolutströmung! — 
2V ——— 
ne ie ne Re 


Die Normalgeschwindigkeit außerhalb der Platte, aber in der Plattenebene wird 


Wera arts) sr 


worin nach (16) A—= R?—.c? zu setzen ist. Ersetzt man noch den arc tg durch den arc sin, so 


erhält man (1): ; > 5 ER 
Du Li en a.) ee FEED): 


Auf der Platte selbst wird v„ natürlich Null. Der Impuls der Flüssigkeit, der zur Berechnung 


der Stoßkraft berücksichtigt werden muß, ist durch das über die eine Seite der Platte erstreckte 
Integral: 


2! 4 
= * m um . — __ 
B o | @-ar ser [ve MıRdRB=- Ve. 2.2...) e 


gegeben, wobei (17) benutzt wurde. Für die Kraft auf den Körper mit dem momentanen 
Stoßkreisradius c folgt damit 


4 d 
ee ea 
—al ) u Se ar 
Dies ist also die Kraft auf den eintauchenden Körper zur Zeit t, wenn er die Eintauchgeschwin- 


digkeit VA) besitzt. Schlägt der Körper zu Beginn des Eintauchvorganges mit der Geschwin- 
digkeit V, auf die Wasseroberfläche auf, so gilt, wenn & seine Eintauchtiefe bedeutet 


> 4 d 5 4 4 5 
mi-=— ?=--30,(N)>m,  a)=zeeayerrr 0. (22) 
dt 5 Ä 
und wegen u(c) = re) Z, auch 2=ui(c):c so daß 
50 V, ER VER 40 
ET - eye m) nl GEeEe. .0 Le ue . (23), 
3m 
„ OR, ° 3 Ugc? 
= — 1, Zu to el 0 SEE 
(Fo ÄRZTE mE le) Sc; 


R- 


A 
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Hieraus folgt sofort die Darstellung von P durch c: 
= 3 WwmV3c? 
we) (1 + Me)? 
Ersetzt man in (23) & durch «(c)c, so erhält man durch Integration die Zeit t als Funktion 
von c: 


ER. (25). 


c 
Be len elegant en siksener (26). 
0 
Da endlich die Eintauchtiefe x durch das Integral: 
c 
«= [ u(e) de nice ib aereern > (27) 
0 


ebenfalls als Funktion von c gegeben ist, liefern die drei letzten Relationen den Zusammenhang 
zwischen Stoßkraft, Eintauchtiefe und Zeit in Parameterform. Ist also für eine gegebene 
Kontur die Größe u(c) berechnet, so ist bis auf Integrationen alles geleistet. 

Zum Schluß mögen diese Formeln für die drei oben behandelten Sonderfälle ausgewertet 
werden. 

1. Der gerade Kreiskegel mit dem halben Öffnungswinkel «a: 


sr _ 12 ee 
u(c) = 7 ctg ao; 12 = WmtgaV? ee 
=. Pe e. I. Ma al 
= 7 etga ER dee leı, 7.0: 
2. Die Kugel vom Radius Eins: 
Bleche Im ©, a 2 äh 
AO, = ma: EN MEN TEEETE A me)‘ 
se Le en. e Ko|t 2 Ar hair © od ) 
= z11 ee Min, sec 3), +3e ce. 
3. Der Sinus-Kopf: 
4K(d)— E(e) _ 3MomVo, e N 
N u NOS 


„2 Bo-U- RO). 


Yt= «+ a) El) — 21— 0)K(o)]. 


Bemerkung zur Fehlerabschätzung beim Adamsschen 
Interpolationsverfahren 
Herrn R. v. Mises zum 70. Geburtstag gewidmet 
Von W. Tollmien in Göttingen 


Im Jahre 1930 veröffentlichte R. v. Mises [1] in dieser Zeitschrift eine Abhandlung 
„Zur numerischen Integration von Differentialgleichungen‘“, die sich für die Entwicklung der 
Praktischen Analysis als bedeutsam und folgenreich erwiesen hat. Einmal enthält dieser 
Aufsatz einen nachdrücklichen Hinweis auf das Adamssche Verfahren zur Integration ge- 
wöhnlicher Differentialgleichungen, das zu jener Zeit in weiten Kreisen noch unbekannt oder 
zum mindesten nicht nach Gebühr gewürdigt war. Das ursprünglich Adamssche Verfahren, 
das auf einer fortlaufenden Extrapolation beruht, wurde in der anschließenden Entwicklung 
verbessert, indem Interpolationen zur Sicherung und Korrektur der Extrapolationen heran- 
gezogen wurden. Auch eine schematisierte ‚‚Anlaufrechnung“ wurde entworfen, um nach 
festen Regeln die Ausgangswerte zu erhalten, die zur Ingangsetzung des Adamsschen Ver- 
fahrens notwendig sind. In jener Abhandlung [1] hatte R. v. Mises des weiteren zum ersten 
Male eine Fehlerabschätzung für das Adamssche Verfahren gegeben, welche die Fehlerfort- 
pflanzung voll berücksichtigt, während früher nur der Fehler beim Einzelschritt betrachtet 
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worden war. 1938 schrieb ich selbst einen kleinen Aufsatz [2], in welchem ich das Problem 
der Fehlerabschätzung von einem anderen Standpunkt aus angriff. Die von mir unter einem 
gewissen Vorbehalt vorgeschlagene Fehlerabschätzung war günstiger als die von R. v. Mises. 
Da bei der Durchführung meiner Fehlerabschätzung noch eine gewisse Schwierigkeit bestehen 
geblieben war, auf die ich gleich zu sprechen kommen werde, setzten nach geraumer Zeit 
Bemühungen von anderer Seite ein, den Abstand zwischen der Misesschen Fehlerschranke 
und der von mir vorgeschlagenen zu verringern. J. Weißinger [3] verschärfte die Misessche 
Fehlerabschätzung, während E. Mohr [4] sich enger an meinen Gedankengang anschloß. 


Der siebzigste Geburtstag von R. v. Mises gibt mir nun einen willkommenen Anlaß, 
auf diese Angelegenheit zurückzukommen und kurz darzulegen, wie ich selbst mir die Behebung 
jener Schwierigkeit vorstelle, die in meiner früheren Arbeit noch enthalten war. Mein Aus- 
gangspunkt war eine sehr gut bekannte Abschätzungsformel aus der Theorie gewöhnlicher 
Differentialgleichungen. Es sei eine Näherung %(&) für das Integral y(x) der Differential- 
gleichung 


da 
= fa) on ee ee ee 
bei der Anfangsbedingung %(2,) =Y, gesucht, wobei die Näherung der Differentialgleichung 
dy \ 
Fe Y) + Bla)r. ae 


genügt. Ist ö eine Schranke für den Fehlerbetrag |y, — %| bei 2 = x,, so gilt 
- M |e—e;. € , M|a—z;| e 
|y — Y|< öde ei| + y(e OT eh ie a OR (3% 


M ist die Lipschitzsche Konstante für f: 


If(®, y) -f® y)l= My —Yıl -» >: 000.0. (A). 
Die Definition von e erfolgt durch die Beziehung: 


Ts) dal <ee—a,| eh DE DE u 
© 


| 


oder etwas roher durch die Beziehung 


Auf die Anlaufrechnung beim Adamsschen Interpolationsverfahren will ich hier nicht ein- 
gehen (vgl. dazu etwa [2]). Die Näherung % werde für äquidistante Argumente , =, +ih 
berechnet, die wir als Stützstellen bezeichnen wollen. A ist die Schrittweite. Es werde 


abgekürzt: 
ya)=h, Sm Ia)=Fa)=F .....:2222.lM. 


Verwendet man ein Interpolationspolynom 3. Grades für F, so erfolgt bei der ‚fortlaufenden 
Rechnung‘ (nach Abschluß der Anlaufrechnung) die Berechnung der Näherungen nach der 
Formel: 


£ x 1 1 en 
Yn+1 = nt ru 78) Vn+ıf — 13 Varnf- 54 Vn+1 on...) 
Die Näherung % genügt also der Differentialgleichung: 
By 
7% = F(x) Tr era Di a er (9), 


wo F(z) das Interpolationspolynom 3. Grades ist, das in den Stützstellen &,41, &, &n-1, 
%,_2 die entsprechenden Werte F,+1, F,, Fn_ı, F„-2 annimmt (s. Bild). Gemäß den For- 
meln (1), (2) und (9) ist 

d(&) = FR ala Ra SE 


also bis auf das Vorzeichen der Interpolationsfehler bei der Approximation von F. Sobald 
man diesen Interpolationsfehler in brauchbarer Weise abgeschätzt hat, so ist die Bestimmung 
des Methodenfehlers bei der fortlaufenden Rechnung des Adamsschen Verfahrens nach (3) 
erledigt, wobei ö die Schranke für den Fehlerbetrag am Schluß der Anlaufrechnung ist, die 
bei & = 2, enden möge. 


2. } th. 1 . > . < 
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Bei der Abschätzung des Interpolationsfehlers [vgl. (10)] ergibt sich nun die erwähnte 
Schwierigkeit. Man würde zunächst geneigt sein, das Cauchysche Restglied für den Inter- 
polationsfehler zu benutzen, also 


De > 2 m mr nn) a) 5 


. (11) 


setzen, wo & ein Zwischenwert zwischen &(2,<&< %,+1) und den 4 Stützstellen &,+1, 


ns Sn—ıs &n-2 ist. Jedoch ist F(x) in den Intervallen zwischen den Stützstellen durch 
verschiedene Polynome 3. Grades definiert, z. B.: 


m UL N ei) P} 1 } 
Fa) Far + #7 ,r+@ > u : ah 
E i . (12a), 
— %n U — in) U — In 3 « a > 
ar Ki +) ' 3! = ( S ) VnrıF für Ln Zeh Lu+1 
dagegen 
7 i2 UT Un (X er, %n) (® == %.-1) 2 (®% FU %,) (2 ae Te) (® ar, Un—2) 3% 
ACH Ra Na ern 21 hr Var+ 311 Vnf 
für ER ee A (12b). 


Kn-3 2 r -] ° A x 


Der Verlauf von F über zwei Schrittweiten ist durch eine aus- 

gezogene Linie angegeben. Die gestrichelten bzw. punktierten 

Linien stellen die beiden verschiedenen Interpolationspolynome 

dar. Die durch Nullenkreise markierten Punkte entsprechen den 
Werten F(@,). 


Daher springt bereits die 1. Ableitung von F(«) an den Stützstellen (s. Bild). Würde man 
nun formal bilden: 

ren), tan, a fer 

N 9209 de 9Y?\de dy da? 


und in diesem Ausdruck nach (9) 2 durch # ersetzen, so sieht man, daß bereits die 2. Ableitung 


von F an den Stützstellen nicht mehr existiert, erst recht natürlich nicht die 4. Ableitung!). 


Ein Ausweg aus dieser Schwierigkeit bietet die Verwendung des Amp£reschen Rest- 
gliedes für den Interpolationsfehler 


— dx) —z F(&)—F(«) = (C— &+1) (® Zu %n) (2— %—ı) (%— 2,2) Fix, Un+1, Ins Un—1» Cn—] 
%n <TF< Dur (13). 


Die 4. dividierte Differenz am Schluß dieser Formel existiert. Dieser Vorschlag ist schon in 
meinem Aufsatz [2] von 1938 enthalten. Leider gibt es bislang keine bequeme Möglichkeit, 
diese 4. Differenz unmittelbar abzuschätzen. Daher möchte ich jetzt ein anderes Vorgehen 
empfehlen. 

Dabei werde ausgenutzt, daß y und F=f(z, y(x)) nur an den Stützstellen tabellarisch 
bekannt sind, daher in den Intervallen zwischen den Stützstellen leicht abgeändert werden 


1) Diesen Sachverhalt hatte ich schon im Sommer-Semester 1939 in einer Vorlesung über ‚Praxis 
der Differentialgleichungen“, die ich an der Technischen Hochschule Dresden hielt, eingehend behandelt, 
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können, wenn man nur dafür sorgt, daß an den Stützstellen die Näherungswerte für y und 
(= y) erhalten bleiben. Wir führen also eine neue Funktion %* ein, für die 


va) ya), Va) Rla), dentimn—-ın-2.....(1) 


gelten soll. Naheliegenderweise nahmen wir für diese neue Funktion y* das Interpolations- 
polynom 7. Grades P,(x), das gerade durch die 8 Bedingungen (14) bestimmt ist. Zum Beispiel 
liefert die Lagrange-Hermitesche Interpolationsformel: 


n+1 5 n+1 R , K 
P= 2 au 2 ra er ak ONE 


i=n—2 


wobei die L, die Lagrangeschen Elementarpolynome 3. Grades bezeichnen, also z.B.: 


(% — %n-2) (®— %) (& — In+1) 

(&n-ı ae N N .; %n) 9,24 u Cn+1) 3 

Statt dessen könnte man auch ein Newtonsches Interpolationspolynom mit zusammen- 
gerückten Stützstellen nehmen, Die so festgelegte Funktion y*—=P, wird ausschließlich für 
das Intervall ,<x&<&,+ı in Gedanken als neue Näherung des Integrals y% von (1) ge- 
nommen, während in dem vorhergehenden Intervall &,-2<x<x, das Polynom y*=P, 
nur herangezogen wird, um eine bequeme Abschätzung des Interpolationsfehlers zu erhalten. 

Weiter werde in dem Intervall 2, <x<s<%,+1ı die jetzt grundsätzlich bekannte 


In ee 


Funktion 
* 

eingeführt. Da - — Pr ist, so ist also jetzt der Interpolationsfehler 
—d=-F*—F, 

für das exemplarische Intervall &,<&< &,+1 abzuschätzen, nachdem analog zu (2) 
dy* 
— =fF*+% 
dx 


für das Intervall ,<x<@,+ı gesetzt ist. In bekannter Weise wird für diesen Inter- 
polationsfehler der Ansatz gemacht: 


— Bd = F* (x) — P;(a) = (8 — 1u—2) (° — 8-1) (E — %) (®— Inr1)X(@) - - (17). 
Zur Abschätzung dieses Fehlers, der, wie gesagt, für das Intervall ,<x<2%,+] zu be- 


stimmen ist, wird nun das größere Intervall &,_2*-- & +1 benutzt, also F* in diesem grö- 
Beren Intervall gleich f(x, P,(x)) gesetzt. Die Hilfsfunktion: 


p(A) = F*(A) —{Pz(A) + (A— m.) (A— n-1) A— m) A— mr) - - - (18) 
verschwindet für A=&,_3, %,—1, Ins %n+ı und für A=z, also an mindestens 5 Stellen. Wegen 
der Stetigkeit der Ableitungen von F ist die viermalige Anwendung des Rolleschen Satzes 
möglich, wodurch folgt, daß gV(}) mindestens eine Nullstelle bei A= xzhat, wo x zwischen 
Dn—2, Cn—1, Cy, &n+ı und liegt (x selbst liegt zwischen x, und &,;1). Da 

FA) = PREV) — PR) — 41 4(a) 
ist, so gilt für A= x: 
0 = FIVE) — PV (zZ) —4!1x(e). 
Also ist 
FW (2) — PV)(z 
en 1 errukteheheg Arch. Mau lo 


wobei &%,-—2< % < %,+1 ist. Führt man in der Faktoriellen (& — # „-.) (2 — %.-ı) (C— %,) 
(2 — %,+1), die in (17) auftritt, &= 2,_12 + uh ein, so wird diese Faktorielle gleich 
h* (u as) Bu ) und ihr Höchstbetrag, der bei = — angenommen wird, ist ht, wie 


4 4 
man leicht feststellt. Als Schlußergebnis erhält man 


; 
8] = Fun) — PN)...» - el 
Das Glied P(V(%) enthält im wesentlichen die Korrektur an einer entsprechenden Formel, 


welche man mit einer bedenkenlosen formalen Anwendung der Cauchyschen Resgliedformel 
erhalten würde. P(V(z) ist ein Polynom 2. Grades. 


Ö 
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Wenn man d gemäß (20) abgeschätzt hat, so kann man die Formel (3) sinngemäß auf y* 
anwenden, womit die Fehlerabschätzung geleistet ist. 6, die absolute Fehlerschranke für die 
Ausgangswerte, die aus der Anlaufrechnung zu bestimmen sind, läßt sich ohne Komplikation 
durch einfache Anwendung der Cauchyschen Restgliedformel herleiten. %k ist bei dem in 
dieser Betrachtung zugrunde gelegten Approximationsgrad gleich 3. 

Es sei betont, daß auf Aussagen über den Fehler einer Näherung innerhalb der Intervalle 
zwischen den Stützstellen verzichtet wird; Feststellungen werden nur getroffen über den 
Fehler der tabulierten Näherungswerte y; an den Stützstellen 2;. 

Die praktische Ausgestaltung dieser hier skizzierten Fehlerabschätzung sowie ihr Ver- 
gleich mit den Ergebnissen von J. Weißinger [3] und E. Mohr [4] sei einer anderen Gelegen- 
heit vorbehalten. 
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Verallgemeinerungen des Seidelschen Iterationsverfahrens”) 
Herrn R. v. Mises zum %0. Geburtstag gewidmet 
Von Johannes Weissinger in Hamburg 


$1. Einleitung 


In [1] haben v.Mises und Pollaczek-Geiringer u.a. gezeigt, daß das lineare 
Gleichungssystem 


n 


AUF TEE A REN. Se serieL) 
k=1 
oder in Matrizenschreibweise kurz 
Ge R . .L.2) 
nach dem Seidelschen Verfahren der Iteration in Einzelschritten, d.h. durch 
1 IS n | : ; 
art) = — — > Yerotl) 4 >. a) — rl s Dun 2 re ll,n) 
a A e! ki+l 


mit beliebigen Anfangswerten a gelöst werden kann, wenn die Matrix Y symmetrisch und 
definit ist. Allgemeiner ist dann auch eine ‚‚Iteration in Gruppen‘ möglich, ‚‚bei der man 
zur Verbesserung der Unbekannten &,, 25, %,, ... die Gleichungen mit den Nummern «, ß, y 
verwendet und mit «, ß, y,... in zyklischer Folge sämtliche Gleichungsnummern durchläuft“. 

In 83 der vorliegenden Arbeit wird zunächst dieser Satz für eine etwas allgemeinere 
Klasse von Iterationsverfahren und Matrizen im Rahmen einer allgemeinen Fixpunkttheorie [3] 
neu hergeleitet. Der benötigte Fixpunktsatz sowie ein grundlegender Hilfssatz über Matrizen 
ist in $2 zu finden. Anschließend wird der Satz auf symmetrische, semidefinite Matrizen 
verallgemeinert. Aus dieser Verallgemeinerung kann man dann leicht die (determinantenfreie) 
Theorie beliebiger linearer Gleichungssysteme gewinnen. 

Die Gedankengänge des $3 lassen sich leicht auf allgemeinere Probleme übertragen. 
In $4 wird die Übertragung für lineare Fredholmsche Integralgleichungen zweiter Art 
durchgeführt. Als Ergebnis erhält man die vollständige (determinantenfreie) Auflösungstheorie 
dieser Integralgleichungen. 


$ 2. Beweishilfsmittel 


Gegeben sei eine additive Gruppe 2 von Elementen £, n,...., für die eine Norm |&] mit 
folgenden Eigenschaften definiert sei: 

a) Es ist |-£]=]£|>0, und |&| = 0 ist äquivalent mit &=0. 

b) Für beliebige &, n gilt |E + nl<I&l+Inl- 


*) In dieser Arbeit werden Teile eines auf der Freiburger Tagung der GaMM im März 1951 gehaltenen 
Vortrags ausführlicher dargestellt (vgl. [4]). 
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c) Die GlSESENE lern &,|—=0 für eine Folge £, ist gleichbedeutend mit der Existenz 
eines Grenzpunktes & — _ Jim &, (d.h. im lE— &,| =). 


vo 


Ferner sei 7 eine additive AR 
T(e+m)=TE+Tn 
für beliebige &,n, die 2 in sich abbildet. Existiert eine reelle Zahl |7]| mit 
ızEeI<]|7|jE| für alle £&, 
so heißt |7T| ein Betrag von 7. Dann gilt der 
Fixpunktsatz: Besitzt die Abbildung T einen Betrag |T|< 1, so hat die Gleichung 
a nie 


bei gegebenem n genau eine Lösung E, und diese ergibt sich als 


E=lim &, 
aus der Folge f 
,auent Ts vl a, 
mit beliebig gewähltem &,. Es gilt die Fehlerabschätzung 
ıT| Kar 
|E ws Ele 1 = 7 | =: HS; ur 2 &,| ee ee (2.3). 


Anmerkung: Gelegentlich hat man einen En. der die Eigenschaften a), b), 
nicht aber c) besitzt, so daß 2 nicht vollständig ist, eine konvergente Folge &, 
(d.h. ae | &,|=0) also nicht unbedingt ein Grenzelement besitzen muß. Dann folgt aus 


7] < 1 "Zunächst nur die Konvergenz der Folge (2.2) sowie die Existenz höchstens einer Lö- 
sung von (2.1). Kann man aber für die Folge (2.2) die Existenz eines Grenzelements nach- 
weisen, so bleibt der Satz in vollem Umfang richtig. 

Hinsichtlich des sehr einfachen Beweises sei z.B. auf [3] verwiesen, wo der Satz in 
allgemeinerer Form formuliert und bewiesen ist. 

Grundlage der folgenden Betrachtungen ist eine Beziehung zwischen quadratischen 
Matrizen, die mit großen deutschen Buchstaben bezeichnet werden sollen. Unter W werde 
die Transponierte von X verstanden. & bezeichnet die Einheitsmatrix. 

Hilfssatz: Ist © eine symmetrische, B, eine beliebige Matrix derart, daß (S— B,)! 
existiert, und setzt man mit beliebigem %: 


B=-+F AI-6-9—-B C=-(5-B)1B:-..: . . (24) 
so galt 
CAE—-A—- (E-C)(S+FH(E-O HB TEE) . . . . . (25). 
Beweis: 
Es. ist 


AE-(S-H—- B)E=-B(E- ©). 
} P=-(E-E)SHFlE—E ..... 20. (26) 
wird daher 


CAT —- — (E-V)B(E- G+ATC- — B-(E-C)(B,— S- FIE-E)+AC 
— — BP+(E-C)(S— B)(E-E+ ACT — PHS- ÜH—- B)(E-O+LAE, 


Mit der Abkürzung 


also 


TAC-— BPHUAU- WE-OHN, 
Voy zen 


woraus wegen 
die Behauptung folgt. 


$ 3. Lineare Gleichungssysteme 
Satz 1: VW und S+% seien positiv definite, symmetrische Matrizen, es sei 
Y=-6- Yı-, 1 BR SSH Fe (3.1) 
und es existiere (S— B,)!. Dann konvergiert die mittels 


=(S—- 8)1%, w=(© Die (3.2) 


Ö 


> 2 
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gebildete Vektorfolge 


Yrı=lm+tr, Ver len ee (3,3) 
bei beliebigem x, gegen die (eindeutig bestimmte) Lösung von 
ei aaa) nee erh ben (BR): 
Beweis: Offenbar ist jede Lösung von Ar =r auch Lösung von 
ae aan. N 


und umgekehrt. Man kann also in der Behauptung auch (3.5) an die Stelle von (3.4) setzen. 

Wir zeigen, daß der Satz ein Spezialfall des Fixpunktsatzes ist. Als additive Gruppe 2 
wählen wir die Menge aller n-dimensionalen Vektoren t mit der gewöhnlichen Vektoraddition. 
Als Norm nehmen wir den Betrag, der sich aus der zu X gehörigen Metrik ergibt und den wir 
zum Unterschied von dem gewöhnlichen Vektorbetrag |r| mit <r> bezeichnen wollen: 


ED aan. rn ea (8.6). 
1 


Bekanntlich ist hierdurch wirklich eine Norm mit den Eigenschaften a, b, c des $2 definiert. 


Als Abbildung &*—=TE nehmen wir die lineare Transformation x*—€r, die offenbar 
additiv im Sinne des $2 ist. Der Satz 1 ist also bewiesen, wenn die Existenz einer Zahl |C]< 1 
nachgewiesen ist, so daß für alle r gilt 


<Er> — VrVAC:< Cl <r> ee (320). 
Aus den Symmetrievoraussetzungen und (2.7) folgt die Symmetrie sowohl von © wie 

von %. Daher ist nach dem Hilfssatz: 
Be Na en... .(38), 


wobei ® nach (2.6) wegen der positiven Definitheit von S+% positiv definit ist. Mithin gibt 
es eine positive Zahl 9, so daß 


BIT mer Deren... 20 2. (3,9) 
gilt. Folglich ist (3.7) mit |C||=Y1—p<1 erfüllbar und der Satz bewiesen. 
Bemerkungen: 


1. Offenbar bleibt Satz 1 richtig, wenn man den Begriff ‚‚positiv definit‘ durch ‚‚negativ 


definit‘“ ersetzt. 
2. Zerlegt man die Matrix X nach folgendem Schema: 


© 


d.h. wählt man für © die Hauptdiagonalmatrix von X bzw. eine aus lauter ‚„‚Hauptminoren“ 
bestehende Matrix, während — ®, (—%,) aus den übrigbleibenden Elementen unter-(ober-) 
halb der Hauptdiagonale besteht, so wird die Vektorfolge x, des Satz 1 identisch mit der 
nach (1.3) gebildeten Folge r, des Einzelschrittverfahrens bzw. mit der entsprechenden Folge 
des Gruppenverfahrens. Satz 1 umfaßt also die in $1 genannten Ergebnisse aus [1]. 

3. Für eine beliebige Matrix X wollen wir in dieser Bemerkung unter * stets die sym- 


metrische Matrix 
Yr—4(A+HW) 


verstehen. Die zu Yund X* gehörigen quadratischen Formen sind dann gleich, so daß Definit- 


heit von UA mit Definitheit von A* äquivalent ist. 
Geht man in der Relation (2.5) des Hilfssatzes links und rechts zu den transponierten 
Matrizen über und mittelt die entstandene Gleichung mit (2.5), so entsteht 


CA = U — Pr — Dr, 
P-=(E-E)(S+HN(E-G ud = HE), 


wobei 


also 
Br = (E-)(S+F)LEC— €) 
IH —- H(E-G) +HE-C)E-HD)=-IWP-WEHEAN—-W) 
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ist. Fordert man demnach in Satz 1 an Stelle der Symmetrie von YWund S-+7% (und damit 
von W, ©, %) lediglich die Symmetrie von ©, so bleibt, wie der Beweisgang sofort erkennen 
läßt, die Behauptung, insbesondere z.B. auch die Konvergenz des Einzelschrittverfahrens, 
bestehen, falls B*+U* positiv definit ist. Aus Stetigkeitsgründen ist diese Voraussetzung 
für nicht allzu unsymmetrische, positiv definite Matrizen W, für die also die Elemente von 
A — X absolut genommen genügend klein sind, sicherlich erfüllbar. Da wegen 
Pr +O9r=(E-O)RE— 6) 
=-6+3 +15 - HITZE SIE ee 

PB* + * gleichzeitig mit NR positiv definit ist, kann man diese Eigenschaft auch von R ver- 
langen. Hinreichend dafür ist |R— E)1<, wobei |W| einen der üblichen Matrizenbeträge 
(vgl. [3]) bedeutet. Ein ähnliches Kriterium ist in [2] hergeleitet worden. 

Nach diesen Bemerkungen wollen wir jetzt den Fall untersuchen, daß WA symmetrisch 
und semidefinit ist, daß also nur Y’ Ur >0 statt Y’ Ar >0 für alle r gilt. Dann sind die Diagonal- 
Elemente a;; von X nur noch > Ostatt >0, so daß das Einzelschrittverfahren in der Form (1.3) 
u. U. gar nicht durchführbar ist. Dagegen existieren stets Zerlegungen von Y, die alle Voraus- 
setzungen des Satz 1 erfüllen. Bezeichnet nämlich ® die Diagonalmatrix von W, B, bzw. 
B,—®B; die links bzw. rechts der Diagonale stehende Matrix: = D+B,+B,, so braucht 


man nur 

s=-d+% ıı=- %=- 845% 
mit irgendeiner symmetrischen, positiv definiten Matrix % zu setzen. Wählt man z.B. =, 
so erhält man statt (1.3) die Iterationsgleichungen: 

[SS 


n 
WAı — — In. ; WM) —_ 9, — F; — ae . . 
%; a; . ET DE + Dan at %; , ı R 2, 0% n (3 10) 


k=i+1l 


Die Vektorfolge £,41 =Ct» +t* ist dann für beliebiges x, definiert, aber sie braucht nicht 
zu konvergieren. Denn im Konvergenzfalle muß sie offenbar gegen eine Lösung von t—= &r +1r* 
d.h. von Ar=tr konvergieren, und bei singulärem W existiert eine solche Lösung gar nicht 
für aller. Wir werden aber sehen, daß das Verfahren immer konvergiert, wenn Ar —=tlösbar ist. 

Wir betrachten eine Zerlegung der symmetrischen, positiv semidefiniten Matrix U wie 
in Satz 1 und führen folgende Bezeichnungen ein: 

R, ist der Raum der Nullösungen von X, besteht also aus allen Vektoren n, die der 
Gleichung An=0 bzw. damit äquivalent n—=&ngenügen. R,ist der zu R, orthogonale Raum, 
und ,=(& — B)"!R, besteht aus allen Vektoren t= (& — B,)-!3 mit SC R,. 

Hilfssatz 1: Der Raum R, wird durch & in sich abgebildet: CR,< R,. Durch 


Cr VIERTE ee ee a 
wird eine Norm in R, definiert. 
Beweis: Wegen 
CR =KS—- B)TR=-(S—-B)TB(6— B)!R, 
ist die Behauptung ER,<R, gleichbedeutend mit 
B(S—-B)TR,ch, 
und dies wieder mit 
(6 -B)'B%r=0 
für alle n< R, 8CR,. Das ist aber wegen Y’n=0, also n = (6 — 3)! Ban richtig. 
Offenbar ist durch (3.11) eine Norm in R, definiert, wenn Xin R, nicht nur semidefinit, 
sondern sogar definit ist, wenn also t’Atin R, nur für t—=0 verschwindet. Da Y Ar =0 äqui- 


valent mit Ar—=0 ist, haben wir also zu zeigen, daß R, und R, nur den Nullvektor gemein 
haben: Rn R,=0. Ist nun 


(SS — B)3=n, (SCH, ncKR,) 
ein solcher gemeinsamer Vektor, so gilt 


0=1r3=-1 (6 — ai 
S+% ei, 
u 3 n+twon=n =; 1 


Da © +7 positiv definit ist, muß also n—0 sein. 


Ö 
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Nunmehr ergibt sich leicht der angekündigte 


Satz 2: Zerlegt man die symmetrvsche, positiv semidefinite Matrix U wie in Satz 1, was 
stets möglich ist, so konvergiert die Folge 


Daı=ör, te, v—=(), Ih, PR TR 


bei beliebigem Yu < R, gegen eine Lösung von Ug—r, wenn diese Gleichung überhaupt Lösungen 
besitzt. 


Beweis: Da rc R, notwendig für die Lösbarkeit von Art ist, liegen X = (&— B,)"!r 
und wegen ER,c R,alle x,in R,. Da YA und damit auch &— Ein R, nur den Nullvektor annul- 
liert und © +%% nach Voraussetzung sogar im ganzen Vektorraum positiv definit ist, sind 
P=(EC-C)(S+F)(E—E) und A in R, positiv definit. Der Beweis von Satz 1 läßt 
sich also wörtlich übertragen, wenn man überall den ganzen Vektorraum durch R, ersetzt. 

Wir wollen jetzt zeigen, daß die Einschränkung t,< R, unnötig ist, und das Verhalten 
der Folge t, für ganz beliebiges x, und rt untersuchen. 


Hilfssatz 2: Die Folge 
w=tyur=l, 2,...) Y beliebig, 
vst stets konvergent. y— limy), liegt in R, und durchläuft ganz R,, wenn y, den ganzen Vektorraum 
vo 


durchläuft. Die drei Aussagen 

a) 53 y»=3&0’y, konvergiert, 

v=0 v0 

b).9=0 und co) <R 
sind äquivalent. 

Beweis: Da 3— X, für beliebiges Y, in R, liegt, so gilt 

ze DIISEZ OU NM EH I ea 12), 

so daß also nach Satz 2 mit t*=t, 1,=0 die Folge 


—1 
= 2t-(E-©)y=y— ho, =]; ZUEER 
k=0 


konvergiert. Daher existiert auch 9=lim y, und genügt wegen Yy+ı =Ch, der Gleichung 


v— co 

y=6h, d.h. Ay=0. Ist umgekehrt n< R, (also n—=En) vorgegeben, so wird lim y,—1, wenn 
o= 1 gewählt wird. vo. 

Wegen ER,< R, liegen alle x, und damit auch lim =, —Y in R,. Die Aussage y,c< R, 


YVr& 
ist daher gleichbedeutend mit 9, < R, und folglich — wegen y< R,und Run R,=0 — mit y—=0. 
Ferner folgt ‚aus ),<R, wie eben nach Satz 1 die Konvergenz von I &’y, = &Y,, und um- 
v=0 v=0 
gekehrt hat die Konvergenz dieser Reihe 4 = lim 4, —0 zur Folge, womit der Hilfssatz voll- 
vy>oo 


ständig bewiesen ist. 
Wir betrachten nun die Folge 


Deal er all Au 

mit beliebigen Vektoren x, und 1r*—=(6&—B,)!t. Durch Differenzenbildung entsteht 

y=ly-ı =%%, v=tırı — 127, 
mit 

yeh-hb=rr— (En. 
Da nach (3.12) der Vektor t=(€—&)r, für beliebiges tr, in R, liegt, so ist {* CR, (d.h. rc R,) 
äquivalent mit ),< R,. Nach Hilfssatz 2 ist also P3 Y, Zi (t»— x) und damit die Folge t, 
dann und nur dann konvergent, wenn v< R,. Liegt rt nicht in R,, so ist 

) 2 (Ev+ı — %») im Ey on err* 

eine nichttriviale Lösung der homogenen Gleichung Yr=0; läßt man reine Basis des ganzen 
Vektorraums durchlaufen, so liefert dieser Grenzprozeß ein Erzeugendensystem von R,. Man 
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kann also z.B. aus (3.10) durch Differenzbildung alle Lösungen der homogenen Gleichung 
erhalten. 
Wir fassen die wesentlichen, bisher gewonnenen Ergebnisse zusammen in dem 
Hauptsatz: Ist A eine symmetrische, positiv semidefinite Matrix, so gibt es stets 


Zerlegungen 

IA=-6—-d Bed r% 
derart, daß (S—B,)! existiert. © symmetrisch und © +% positiv definit ist. Bildet man aus 
beliebigen Vektoren X, und x mittels 

L=-(5— BR; r= (6 — 9) 7 h 
die Folge 

Cy+1 —=(r, ats N 1, 24 

so konvergiert die Folge yy = I» +ı— I» stets gegen eine Nullösung von U, und zwar dann und 
nur dann gegen eine nichttriviale N ullösung, wenn x nicht orthogonal zu allen Nullösungen von U 
ist. Durch geeignete Wahl von x kann man so Jede Nullösung erhalten. Exwistiert sogar a, Ip, 


so genügt 3 der Gleichung Uz—=r und umgekehrt gilt trivialerweise für jede Lösung ; dieser 
Gleichung: 3=lim %,, wenn %, —3 gewählt wird. Notwendig und hinreichend für die Konvergenz 


vo © 
der Folge x, (und damit für die Lösbarkeit von Ag —t) ist die Orthogonalität von t zu allen N ull- 
lösungen von U. 
Wir betrachten jetzt das Gleichungssystem 


Me=m 


mit einer beliebigen rechteckigen Matrix M. Jede Lösung dieser Gleichung ist auch Lösung 
der Gleichung Ar =t, wo 
AE-UM, = Mn 
as ist. ist symmetrisch und positiv semidefinit und besitzt dieselben Nullösungen wie 
: RM) = RM). Daher gilt auch für die entsprechenden Orthogonalräume: RAU) — R,(M). 
Für die Lösbarkeit von Yr=r ist rc R,(W) notwendig und hinreichend; das ist für 
r— Mm stets der Fall, da MWmcR,(M) für beliebiges m gilt. Die Gleichung Ar =t ist also 
nach dem Hauptsatz stets lösbar. Für die Lösung u kann man aber nicht Mu=m folgern, 
sondern nur Mu—=m-+p mit vCR,(M). Da MucR,(M) für jeden Vektor u gilt, so ist 
mc R,(M)äquivalent mit v—=0. Ist also die für die Lösbarkeit von Mr=m notwendige 
Bedingung m< R,(M) erfüllt, so sind die nach dem Hauptsatz sicher vorhandenen Lösungen 
von Ar =r auch Lösungen von Mr = m. 


$4. Lineare Fredholmsche Integralgleiehungen zweiter Art 
Die Übertragung der vorstehenden Überlegungen auf die Integralgleichung: 


2(0)— | Kior)ek)de=r(o) cl re Fell 


in der wir der Einfachheit halber die gegebenen Funktionen X (o,r) und r(o) als stetig voraus- 
setzen wollen, macht keine Schwierigkeiten, da die meisten Schlüsse allgemein für lineare 
Operatoren gültig sind. 

Um die Analogien deutlicher hervortreten zu lassen, betrachten wir die Funktionen &(o), 
y(o), r(o), ... als Elemente 2,4, t,... im Raume der stetigen Funktionen, definieren 


1 
cy=yr = (0) y(o) do 
und schreiben 
y=sI Try 


an Stelle von 
1 


y(o) = [ K(o,t) x(t) dt bzw. y(o) = [k (7,0) &(r) dr. 


Damit gilt dann z. B. 
sr ne [Ron &(T) y(o) do dr . 


Wird der Einheitsoperator mit & bezeichnet, so kann man also statt (4.1) N 


Yr=tr, USE—R sl BB een 


» 
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In $3 haben wir ziemlich allgemeine Zerlegungen der Matrix X betrachtet, um auch für 
semidefinite Matrizen das Verfahren stets durchführen zu können. Jetzt dagegen wollen wir 
nur die dem eigentlichen Einzelschrittverfahren nach (1.3) entsprechende Zerlegung betrachten, 
die bei Integralgleichungen immer zum Ziel führt. An diesem Beispie] wird dann auch das 
Verfahren für allgemeinere Zerlegungen klar. 

Wir definieren die Kerne 


Das , Be ra 

rg a FIN T>o} 
mit den zugehörigen en Dı bzw. ®,, sodaß also überall X= B, + B,, d.h. £— B+B,, 
gilt. Der Operator EC—%;, der einer Volterraschen Integralgleichung entspricht, besitzt 
ein Inverses: 


Dior) . (43) 


(E-B)I=-6+8, 2=-!%. 
v=1 
Dem Operator 2 entspricht als Kern die Neumannsche Reihe: 
L(e, N  Bir(an) . ae 1251 (4.A), 


in der B® die »-te Iterierte von B, ı a ist (4.2) äquivalent mit 
(E-6&)r=t*, wobei C=-(E-B)'B,, 1 =(E—-B)'r.... (45), 
d.h. jede Lösung von (4.1) ist auch Lösung von 


&(0) — [t0) area) neh (4,6) 

mit : 
C(o,r) = B,(o,r) + / Lloso, Boa lao 2... (AT) 
7*(0)=r(0) + La) 1 (DIR NR EI NE (4,8) 


und umgekehrt. Wenn (4.6) durch eine Neumannsche Reihe gelöst werden kann, so ist 
damit (4.1) durch zweimalige Bildung einer Neumannschen Reihe gelöst. 

Wir setzen den Kern K(o,r) für das folgende als symmetrisch voraus. Der Beweis des 
Hilfssatzes aus $2überträgt sich wörtlich auf Integraloperatoren, wobei wir bei unserer speziellen 
Zerlegung hier und im folgenden noch S=€ und $=0 setzen dürfen. Insbesondere ist dann 
die in $3 stets auftretende Forderung, daß © +7 symmetrisch und positiv definit sein soll, 
sicher erfüllt. Die Forderung, daß U positiv definit sei, besagt für den Kern K(o,r), daß 


II Kto) 200) 200) 00 dr < fa*(0)do 


für alle stetigen Funktionen x(o) gilt; bei Semidefinitheit ist das —-Zeichen zugelassen. 
Ist U positiv definit, so existieren bekanntlich zwei positive Konstanten c,, c, so daß 


Ba Or... = (4,9) 


- für alle stetigen Funktionen x&(o) gilt. Daraus folgt sofort, daß die Beziehung ID bo Yy=0 


für irgendeine Folge r, gleichbedeutend ist mit an Utz U. d.h: mit Konvergenz im Mittel 


(gegen 0) für die Folge &, (co). Da aus Konvergenz im 1 Mittel i im allgemeinen nicht gleichmäßige 
Konvergenz folgt, ist der Raum 2 der stetigen Funktionen mit der Normdefinition: 


Be. Neon... 2. (4.10) 
nicht vollständig, wohl aber offensichtlich mit der Definition 
|x| = max |x(o)| Br: 2.2 (411) 
0<o<1 


Die Übertragung des Beweises von Satz 1 aus $3 liefert also zunächst nur die Konvergenz 


der Folge 
vr aß; to beliebig, wur. here rlrl2) 


im Sinne der Normdefinition (4.10). Nun gilt aber für eine quellenmäßig dargestellte Funktion 


y(o) _ fx(6,r) z(r)dr, d.h.y=Cr 
12 
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nach der Cauchy-Schwarzschen Ungleichung 


voice) | x: (T) dt, = \ max [0:(o, T)dr, 


0<o<1H 
d.h. 
nI<errr, 
also nach (4.9) 
hI<a<rD, =. 
6 


Für die Differenz %,41 — Lu+1ı = (rt, — T,) der Folge (4.12) gilt daher 
|iv+ı — La+ıl S@&tp — Ip); 


so daß also für die Iterationsfolge (4.12) aus Konvergenz im Sinne von (4.10) auch Konvergenz 
nach (4.11), d.h. gleichmäßige Konvergenz und damit die Existenz einer stetigen Grenzfunktion 
folgt. Der Fixpunktsatz ist also auch mit der Normdefinition (4.10) gültig und der Beweis 
des Satz 1 vollständig übertragbar. 

Um die weiteren Ausführungen des $3 zu übertragen, benötigt man einige Eigenschaften 
des Raums 2 der stetigen Funktionen, die im Fall des n-dimensionalen Vektorraums selbst- 
verständlich waren. 

2, sei ein linearer Teilraum von 2, der jedes Element aus 2 enthält, das orthogonal 
zu allen Elementen des zu Q,orthogonalen Raumes Q,ist. Beispiele für einen solchen’ Raum 2, 
sind der Raum R, der Nullösungen von NW, der zu R, orthogonale Raum AR, und der Raum 
R=(6— B)TR,. 

Ist x, eine Folge aus Q,, die im Mittel gegen ein Element x aus Q konvergiert: 
im, — (mw —r)=0, so muß r in Q, liegen. Denn für beliebiges y<Q2 gilt dann 
yo» 

lim H’r,=Y'z, also Yr.= 0 Tlür alleuxr ©: 
yo 

Wir setzen jetzt X nur als positiv semidefinit voraus und betrachten eine beliebige Folge 
T, ca, m 

im y Ay, =0,; u use (keit 


vy>o 


Nach der Cauchy-Schwarzschen Ungleichung ist dann auch aim 1 Eu A dr —0 für jede Folge dv 


mit beschränktem yU3,=<3,3?. Wegen %,1,—=1 ist für ee "Folge 2 — Mr, sogar |3,| 
beschränkt, so daß also wegen der Vertauschbarkeit von Y=&— ft mit $& sicher gilt 

lin, Ayn0.5 Mlyee Sta Ar : 3 «anal; 14), 

v>» 
Da die 4, beschränkt und gleichgradig stetig sind, kann man aus ihnen eine gleichmäßig 
konvergente Teilfolge auswählen; » möge fortan nur die Indizes dieser Teilfolge durchlaufen. 
Dann ist also 

y —nliben, Yy 
vo 


im Sinne des zu (4.11) gehörigen Konvergenzbegriffes, d.h. im Sinne gleichmäßiger Konvergenz 
für die zugeordneten Funktionen. Aus (4.14) folgt nun y Ay —=0 und daraus 


Ay=0,..d.h. pen 0000 Dee 7a 
Ferner ist nach (4.13), (4.15): 
1 u Huah= ao tudy wu %H = ws u nd») = 
und daher 
im iz —y/(w-4)=0, 

d.h. die t, konvergieren im Mittel gegen 4, so daß y in 2, liegen muß. 

Ist also Q, ein Raum, der mit dem Raum R, der Nullösungen von U nur die 0 gemein 
hat, so hat sich ein Widerspruch ergeben: Eine Folge x, < 2, mit der Eigenschaft (4.13) kann 
gar nicht existieren, d.h. es existieren Konstanten c,>0 und (trivialerweise) c, >0 mit 


arr<rAr <art für alle oo 


Auf Q,ist daher Konvergenz nach (4.10) äquivalent mit Konvergenz im Mittel, für Iteralions- 
folgen also sogar mit gleichmäßiger Konvergenz. 
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Man bestätigt nun leicht, daß die Ausführungen des $3 auch für die hier betrachteten 
Integraloperatoren richtig sind. Wir wollen das Schlußergebnis noch in die Sprache der 
Integralgleichungen übertragen: ». 

Bildet man aus dem (evtl. unsymmetrischen) Kern H(o,r) der Integralgleichung 


OEL H(o,r)2() BOHRER rear (4,16) 
den symmetrischen Kern i 
K(o,r) = H(o,r) + H(t,o) — [H(e.o) H(o,r) do 
und die Funktion - 
Ko) = mio) | Hleo)m(e)de, 


so ist jede Lösung von (4.16) auch Lösung von (4.1) und der damit äquivalenten Integral- 
gleichung (4.6), deren Lösungen durch Iteration (Neumannsche Reihe) bestimmt werden 
können. Umgekehrt ist jede Lösung dieser (bei beliebigen m(o) lösbaren) Integralgleichung 
auch Lösung von (4.16), falls diese Integralgleichung überhaupt Lösungen besitzen kann, falls 
also m(c) orthogonal zu allen Nullösungen der zu (4.16) adjungierten Integralgleichung ist. 
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